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Abstract 

The  number  of  countries  that  are  gaining  access  to  space  is  increasing  and  the  drastic 
increase  in  both  computing  and  sensor  technology  has  allowed  the  development  of  micro  (10- 
100kg),  nano  (1-10  kg),  pico  (0.1-1  kg)  satellites  to  be  developed  and  put  into  orbit  to  satisfy  a 
particular  operational,  experimental  or  functional  need.  The  attractiveness  of  the  micro-Pulsed 
Plasma  Thruster  (pPPT)  is  that  the  relatively  simple  construction,  reliability,  specific  Impulse 
(Isp)  of  about  800- 1000s  and  require  lower  power  levels  than  other  types  of  electric  thrusters  like 
ion,  and  hall  thrusters.  These  types  of  thrusters  also  supply  very  small  impulse  bits  on  the  order 
of  10-100pN-s  per  thruster  tube  and  are  better  equipped  for  fine  attitude  and  control  of  the 
spacecraft,  which  may  be  desirable  for  surveillance  or  reconnaissance  purposes.  These  thrusters 
to  have  a  couple  of  drawbacks  to  them  as  well,  they  are  typically  not  very  efficient  at  converting 
the  electrical  energy  into  thrust  with  efficiencies  typically  around  10%.  The  thrust  can  also 
decrease  as  the  pPPT  approaches  end  of  life. 

The  propellant  is  polytetraflouroethylene  (PTFE),  also  known  as  Teflon™,  which  is 
contained  entirely  within  the  thruster,  reducing  the  complexity  of  the  device,  eliminating  the 
need  for  propellant  tanks  or  feed  lines.  The  thruster  operates  by  producing  an  arc  that  ablates  the 
surface  of  the  material,  and  that  material  is  then  weakly  ionized  by  the  arc  and  accelerated  by 
electromagnetic  forces  induced  by  the  arc  producing  thrust.  The  Air  force  Research  Laboratory 
Propulsion  Directorate,  located  at  Edwards  AFB,  CA,  developed  pPPT  initially,  as  a  smaller 
version  of  the  widely  tested  Pulsed  Plasma  thruster.  The  pPPT  has  not  been  tested  as  extensively 


as  the  PPT. 


This  research  is  intended  to  characterize  an  assembly  of  nine  pPPT  thrusters  arranged  in 
a  'Gatling  Gun'  (GG)  pattern  complete  with  drive  electronics  and  connections.  In  doing  so,  this 
research  will  investigate  the  effects  of  firing  pattern  on  performance  of  the  GG-pPPT  in  relation 
to  propellant  utilization  of  each  thruster,  thrust  duration/pattem,  the  plume  discharge  and 
frequency  response  of  the  driving  electronics. 
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USING  A  GATLING-GUN  CONFIGURED  MICRO  PULSED  PLASMA  THRUSTER  AS 


A  MEANS  TO  CONTROL  MICRO  SATELLITES  WITH  EXTREME  PRECISION 

I.  Introduction 

This  research  focused  on  characterizing  and  determination  of  the  optimum  operating 
conditions  of  this  thruster  module.  The  ability  to  operate  these  types  of  thrusters  without 
excessive  carbon  build  and  failure  has  been  an  issue  for  previous  research  done  at  AFIT  on  these 
types  of  thrusters.  To  determine  the  best  operating  conditions  we  need  to  characterize  the 
performance  of  the  Gatling-Gun  micro  pulsed  plasma  thruster  (GG-pPPT)  under  various 
conditions.  However,  it  is  important  to  understand  the  various  types  of  electric  propulsion  (EP) 
to  gain  an  understanding  of  the  similarities  and  differences  between  micro  pulsed  plasma 
thrusters  and  other  EP  systems. 

Chapter  1  discusses  the  goals  and  objectives  of  this  research  and  then  explores  the 
different  types  of  EP  systems.  Chapter  2  focuses  on  the  historical  research  and  results  of  the 
pulsed  plasma  thruster  in  general  and  the  micro  pulsed  plasma  thruster  developed  afterwards. 
Chapter  3  describes  the  equipment  and  methodology  behind  the  research  conducted.  Chapter  4 
presents  the  analysis  and  results  of  the  research,  and  Chapter  5  draw  conclusions  from  the  results 
and  address  future  research  areas  for  the  GG-pPPT  module. 

The  first  objective  is  to  characterize  performance  of  the  GG-pPPT.  The  GG-pPPT  is  a 
new  device  using  nine  individual  pPPTs  driven  by  a  single  set  electronics  package.  Previous 
research  by  Keidar,  Boyd,  Antonsen,  Burton,  and  Spanjers  [1]  indicates  too  high  or  too  low 
discharge  energy  can  lead  to  latent-time  vaporization,  accumulation  of  carbon  deposits  and  non- 
uniform  arc  formation.  This  phenomenon  can  result  in  lost  efficiency  because  as  the  carbon 
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builds  up  more  power  is  required  to  ionize  and  the  cycle  continues  until  the  carbon  shorts  the 
electrodes  preventing  further  operation.  A  secondary  objective  will  be  to  see  how  the 
performance  changes  with  regard  to  thrust,  and  efficiency  by  adjusting  the  parameters  of  the 
single  set  of  drive  electronics.  To  determine  the  characterization  we  will  need  to  collect  thrust 
data,  firing  frequency  of  the  discharges,  firing  pattern,  and  power  level. 

Research  Objectives 

This  research  has  several  objectives.  First,  characterize  the  GG-pPPT  firing  pattern, 
input  power,  and  frequency.  This  objective  will  be  accomplished  through  variation  of  pPPTs 
receiving  electrical  power,  increase  or  decrease  of  electrical  power  supplied  to  pPPTs,  and 
discharge  time  of  the  main  capacitors  to  the  pPPT.  Second,  we  want  to  identify  the  best 
conditions  to  operate,  focusing  on  optimized  efficiency  and  longer  life  of  the  individual  pPPTs. 
This  objective  requires  determination  of  the  operational  envelope  of  thruster.  This  task  will  be 
accomplished  through  graphical  plots  created  with  power,  frequency,  and  discharge  energy,  all 
unique  to  this  thruster,  and  efficiencies  calculated  based  on  the  ratio  of  jet  power  to  discharge 
power.  Third,  produce  thrust  tables  designed  for  more  precise  control  of  a  satellite  through 
calculation  or  collection  of  torque  data. 

Background 

The  desire  to  have  a  satellite  point  precisely  in  the  same  direction  or  change  its  attitude 
by  very  small  amounts  makes  the  pPPT  the  right  thruster  for  the  job  concerning  Intelligence, 
Surveillance,  and  Reconnaissance  (ISR)  applications  or  for  very  precise  north-south  station 
keeping  (NSSK).  The  thruster  also  finds  applications  for  precision  space  experiments  and 
astronomy.  For  instance,  the  detection  of  gravitation  waves  would  require  very  precise 
positioning  relative  to  the  other  satellites  in  the  constellation  to  detect  these  waves  or  ripples  in 
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space-time.  The  waves  are  generated  by  massive  rotating  objects  like  neutron  stars  or  black 
holes  predicted  by  Einstein’s  theory  of  general  relativity  [2],  They  also  find  application  for 
many  other  small  satellites  like  ‘CubeSats’  where  they  could  be  the  propulsion  source  for  attitude 
adjustment  and  alignment.  Electric  thrusters  have  gained  more  use  in  recent  times,  because  the 
manufacturing  technology  for  sensors  and  equipment  has  advanced  to  the  point  where 
components  are  smaller  in  physical  size.  This  propulsion  method  is  attractive  for  CubeSats, 
because  their  small  size  does  not  require  massive  impulse  bits  from  chemical  rockets  for  re¬ 
orient  themselves.  This  niche  application  for  electric  thruster  propulsion  becomes  more 
attractive  as  precise  attitude  and  control  are  required  for  these  small  satellites  and  sensor 
technology  power  and  size  decrease  with  better  advances  in  materials  and  manufacturing 
engineering. 

The  main  attractiveness  of  the  micro-pulsed  plasma  thruster  is  that  they  are  relatively 
simple  to  construct,  very  reliable,  have  specific  Impulse  (Isp)  of  800- 1000s  and  require  lower 
power  levels  than  other  types  of  electric  thrusters  like  ion,  and  hall  thrusters.  These  types  of 
thrusters  also  supply  very  small  impulse  bits  on  the  order  of  10-200pN-s  per  thruster  tube  and 
are  better  equipped  for  fine  attitude  and  control  of  the  spacecraft,  which  may  be  desirable  for  ISR 
purposes  or  interferometer  requirements. 

They  could  also  find  applications  on  satellites  designed  for  extra-solar  planet  detection. 
Early  detection  methods  used  the  wobble  of  the  star,  and  the  resulting  Doppler  shift  produced  as 
the  planet  orbits  the  star  and  the  system  oscillates  about  its  center  of  mass.  However,  newer 
methods  use  the  dimming  of  start  light  produced  occurring  when  the  planet  transits  the  line  of 
sight  to  the  star,  which  requires  precise  tracking  and  orientation  capability  to  detect.  They  could 
also  counteract  potential  drag  losses  for  satellites  in  Low  Earth  Orbit  (LEO).  The  GG-pPPT  is 
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differentiated  from  other  pPPT  systems  in  that  it  contains  nine  thrusters  and  electronics  all 
housed  in  one  unit,  fitting  into  a  1U  (1000  cubic  centimeters)  package  and  providing  a  smoother 
thrust  profile  than  current  pPPT  devices  in  use  onboard  satellites. 

Electric  Propulsion 

Chemical  rockets  and  thrusters  produce  thrust  through  the  heating,  acceleration,  and 
expulsion  of  the  gas  or  solid  through  a  nozzle.  The  specific  impulse  can  be  thought  of  as  the  fuel 
efficiency  of  the  thruster,  and  for  chemical  rockets,  this  is  limited  to  about  500  seconds  for  a 
Hydrogen-Oxygen  rocket.  Longer  space  missions  require  fuel  efficient  rockets.  This  is  where 
electric  propulsion  has  the  advantage.  Compared  to  chemical  rockets,  electric  propulsion 
provide  a  higher  Isp  resulting  in  a  lower  mass  fraction  for  the  same  change  in  velocity.  The 
other  drawback  to  chemical  systems  is  that  they  rely  on  pressure  forces  caused  by  the  chemical 
reaction  and  process  of  combustion  to  produce  thrust.  Because  the  energy  transferred  to  the 
exhaust  products  comes  from  combustion  of  reactants,  chemical  rockets  are  energy  density 
limited  due  to  enthalpy  limitations,  flame  temperatures  and  temperature  limits  of  the  materials 
exposed  to  these  hot  gases.  Thrust  levels  are  very  high  for  these  devices  on  the  scale  of  a  couple 
of  Newtons  to  Mega-Newtons  (MN). 

Electric  thrusters  get  around  these  limitations  by  using  electro-static,  electromagnetic  or 
electrodynamic  methods  derived  from  the  Lorentz  forces  [3]  and  Maxwell’s  Field  Equations  [4], 
Essentially,  three  classes  of  electric  thrusters  are  employed:  electrothermal,  electrostatic,  and 
electromagnetic.  The  three  categories  produce  thrust  differently,  but  the  main  common  thread 
between  them  is  the  energy  to  produce  thrust  originates  from  electrical  power,  not  the  internal 
energy  of  chemical  reactions  as  the  case  with  solid  or  liquid  thrusters.  We  will  now  examine  the 
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simplest  and  first  kind  of  EP,  the  electrothermal  thruster  that  uses  electricity  to  heat  the 
propellant  and  expand  it  through  a  nozzle  producing  thrust. 

Electrothermal  Thrusters 

The  first  category  of  electric  thrusters  is  the  electrothermal  type  like  resisto-jets  and  arc- 
jets,  which  use  electrical  energy  to  heat  a  gaseous  propellant.  These  thrusters  operate  in  the 
following  way.  First,  the  cathode  is  located  between  prior  to  the  constricted  region  of  the  arc  jet 
thruster  and  this  cathode  produces  a  very  high-energy  arc  that  attaches  to  the  anode  just  on  the 
other  side  of  the  constrictor.  The  equivalent  temperatures  at  the  core  of  these  arcs  are  in  the 
range  of  10,000-40,000  K,  and  typical  wall  temperatures  of  1, 000-3, 000K  [3],  Second, 
propellant  is  passes  along  the  cavity  over  the  cathode  and  through  the  constrictor. 

These  function  very  similarly  to  chemical  rockets  and  have  a  nozzle  that  converts  the 
pressure  forces  to  kinetic  energy  through  expansion  on  the  thruster,  propelling  it  forward. 
However,  the  kinetic  energy  relies  on  the  electrical  energy  to  heat  the  propellant  vs.  the  internal 
chemical  kinetics  of  combustion.  They  are  also  different  from  other  electric  thrusters,  in  that  the 
thrust  and  acceleration  is  a  purely  aerodynamic  momentum  transfer  between  the  accelerating  gas 
and  the  nozzle.  These  thrusters  are  used  on  mini  and  some  micro  satellites  like  Telstar™  and 
have  typical  Isp  of  500- 1000s  with  the  upper  end  requiring  use  of  hydrogen/oxygen  propellants 
[3],  Figure  1  shows  a  schematic  of  an  arcjet  and  resistojet  thruster  [5],  These  thrusters  can 
achieve  better  performance  than  chemical  systems,  but  Isp  is  limited  by  the  temperature 
limitations  of  the  material  while  exposed  to  the  high  temperature  gas,  arc  discharges  and  erosion 
of  the  surfaces  of  the  chamber.  Electrostatic  thrusters  are  the  next  major  class  of  EP  thrusters 
that  produce  thrust  through  electrostatic  repulsive  forces  between  ions  and  an  electric  field. 
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(a)  Arc  jet  schematic. 

Propellant  is  heated  as  it  passes  through  arc  and  accelerated  out  of  nozzle  by  pressure  forces 


Nozzle 


(b)  Resisto-jet  schematic. 

Propellant  is  heated  as  it  passes  resistive  elements  and  accelerated  out  of  the  thruster  by  pressure 

forces  as  it  passes  through  the  nozzle 


Figure  1.  Electrothermal  Thrusters 
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Electrostatic  Thrusters 


These  types  of  thrusters  may  operate  using  electrostatic  forces  to  accelerate  the  particles 
directly  without  relying  on  pressure  forces  to  generate  thrust  as  in  the  case  for  electro-thermal 
propulsive  systems  (resisto-jet/arc-jet).  Ion  thrusters  are  a  common  type  of  thruster  in  this 
category.  The  size  and  complexity  of  these  types  of  thrusters  can  vary  and  power  levels  can 
range  from  a  few  watts  up  to  30  kW  of  electrical  power  for  operation.  Ion  thrusters  work  by 
using  cathode  to  produce  a  stream  of  electrons,  the  electrons  then  strike  propellant  atoms 
(usually  Xe  or  Kr)  injected  through  the  a  hollow  cathode,  or  through  another  port.  The  collision 
of  the  electrons  with  the  neutral  atoms  causes  them  to  lose  an  electron  and  become  positively 
ionized.  An  electrical  grid  accelerate  the  ions  through  electrostatic  forces  with  a  potential 
voltage  difference  of  about  l-2kV  where  they  pass  through  a  grid  and  are  neutralized  on  the 
outside  by  an  electron  beam  to  prevent  charge  accumulation  on  the  thruster  and  hence  the 
spacecraft. 

The  propellants  are  in  the  gas  state  for  these  thrusters  and  has  a  pressure  vessel 
containing  the  propellant  and  are  pressure  feed  systems.  These  thrusters  are  typically  the  highest 
power  consumption  type  of  thruster  due  to  the  high  amount  of  ionization  required  but  offer  high 
very  high  specific  impulse  in  the  range  of  2000-10,000  s  and  efficiencies  approaching  90%  [4], 
Another  type  of  thruster  in  this  category  is  the  Hall  thruster,  and  uses  the  Hall-effect  in 
conjunction  with  electrostatic  forces  to  accelerate  ions  to  very  high  speed  producing  thrust.  This 
type  of  thruster  operates  by  using  a  magnetic  field  or  B  field  that  is  perpendicular  to  the  electric 
or  E  field.  The  magnetic  field  ‘holds’  back  the  electrons,  while  ionization  takes  place,  though 
the  electrons  do  have  a  drift  velocity  and  eventually  reach  the  cathode,  causing  cathode  erosion. 
The  electrons  then  collide  with  the  propellant  ionizing  them.  The  E  field  then  accelerates  the 
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ions  that  have  a  net  positive  charge  out  of  the  thruster,  producing  a  momentum  exchange  and 
producing  thrust  as  shown  in  Figure  2. 

Hall  thrusters  are  electrically  simpler  than  their  electrostatic  counterparts,  but  because  of 
about  30  percent  of  the  discharge  current  produces  electrons,  which  do  not  directly  contribute  to 
thrust,  they  are  about  70%  thrust  efficient  with  regard  to  utilization.  However,  they  are  about  80- 
90  percent  efficient  at  ionizing  the  propellant  so  overall  they  have  a  total  efficiency  between  54 
and  63%  taking  the  products  of  the  two  previous  numbers. 

The  application  of  electrostatic  thrusters  in  deep  space  missions  provides  a  way  to 
explore  our  solar  system  beyond  our  planet  and  the  DEEP  SPACE- 1  spacecraft  is  a  good 
example  of  such  a  spacecraft  [6],  designed  to  test  the  performance  of  an  ion  thruster.  Future 
mission  to  asteroids  could  also  utilize  electric  propulsion  due  to  their  weak  gravitational  field  and 
small  size.  The  last  major  category  of  electric  propulsive  thrusters  is  electromagnetic.  The 
thrusters  use  a  combination  of  electric  and  magnetic  body  forces  to  accelerate  particles  to  high 
velocities,  creating  a  momentum  exchange,  and  producing  thrust. 
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Figure  2.  Hall  thruster  with  external  cathodes,  based  on  Hofer  [7] 

Electromagnetic 

The  electromagnetic  thruster  types  include  the  PPT,  pPPT  and  magneto-plasma  dynamic 
(MPD)  thrusters.  These  thrusters  operate  through  electrical  ionization  of  the  propellant  and 
acceleration  through  external  or  internally  produced  electromagnetic  fields.  In  the  case  of  the 
MPD  thruster,  they  are  really  an  arc-jet  thruster  with  an  electromagnet  surrounding  the  nozzle  to 
aid  in  acceleration  of  the  plasma  and  exhaust  products  through  the  Lorentz  force  once  they 
become  ionized.  A  schematic  of  the  operation  MPD  and  PPT  is  shown  in  Figure  3. 
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Figure  3.  MPD  thruster  operation 


For  the  PPT  and  pPPT,  the  cathode  and  anode  create  a  strong  electrical  field  eventually 
resulting  in  arc  formation  once  the  potential  is  above  the  breakdown  voltage  of  the  material 
(Teflon™).  This  arc  ablates  some  of  the  Teflon™,  surface.  These  particles  pass  through  the  arc 
and  are  ionized.  Self- induced  electromagnetic  fields  accelerate  the  ionized  Teflon™,  electrons 
and  bits  of  cathode  and  anode  out  of  the  thruster.  Through  the  Lorentz  force,  the  thruster 
produces  a  momentum  exchange  producing  thrust  by  pushing  back  on  the  electromagnetic  field. 
The  electrons  and  ions  will  have  opposite  charge  and  will  spin  in  opposite  directions  after 
leaving  the  thruster  but  then  come  together  once  further  away  to  form  quasi-neutral  plasma,  thus 
eliminating  the  need  for  an  external  cathode  to  neutralize  the  exhaust  plume,  as  required  for  Ion 
and  Hall  thrusters.  A  simplified  diagram  of  a  PPT  is  shown  Figure  4. 
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Thrust 


Figure  4.  Schematic  of  Pulse  Plasma  Thruster,  based  on  Wie  [9] 

PPTs  and  pPPT  vary  significantly  in  their  structure  because  of  the  different  sizes  and 

applications.  The  PPT  normally  has  a  spring  fed  propellant  system,  keeping  the  propellant  near 
the  location  of  the  arc  so  it  can  be  ablated  and  ionized  by  the  arc  formation.  The  larger  PPT  also 
has  a  spark  igniter  or  spark  plug  to  reduce  the  amount  of  power  needed  for  initial  arc  creation. 
Once  this  arc  is  established,  it  can  be  maintained  at  a  much  lower  voltage  level  than  the  initial 
trigger  spark.  The  larger  size  of  the  PPT  allows  them  to  be  slightly  more  complex,  but  also 
require  higher  power  levels  than  pPPTs.  PPTs  were  developed  first,  then  pPPTs  afterwards  by 
AFRL  Propulsion  Directorate  at  Edwards  AFB,  CA  in  1997  [8]. 

The  |iPPT  design  is  much  simpler  than  the  PPT.  First,  it  does  not  require  a  spring  to  keep 
the  propellant  in  contact  with  the  arc  face,  because  as  the  electrodes  erode  they  maintain  a  spark 
level  with  the  remaining  electrodes  and  propellant  surface.  The  electrodes  and  propellants  are 
circular  in  cross-section  for  the  pPPT  and  come  in  two  or  three  electrode  variants.  The  two- 
electrode  variant  shown  in  Figure  5  has  just  the  cathode  and  anode.  These  have  since  evolved  to 
the  three  (cathode,  inner  electrode,  anode)  electrode  configuration  for  efficiency  and  power 
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generation  requirements  for  arc  formation.  The  inner  electrode  acts  as  an  improvised  spark  plug 
providing  a  seed  arc  and  plasma  thereby  reducing  the  amount  of  voltage  to  start  arc  creation. 
Once  created,  the  flow  of  electrons  can  be  sustained  at  a  lower  operating  voltage  during  thrust 
generation. 


Figure  5.  Basic  2-electrode  micro-PPT  operation,  based  on  Burton  [10] 

The  relatively  simple  design  and  low  power  consumption  relative  to  other  types  of 
electric  thrusters  makes  them  an  ideal  candidate  for  use  on  some  micro,  nano  and  pico  satellites. 
For  the  PPTs  and  pPPT,  one  of  the  biggest  losses  in  efficiency  comes  from  the  conversion  of 
electrical  energy  into  thrust.  The  ablation  of  Teflon  through  heating  of  the  surface  does  not 
produce  a  strongly  ionized  gas  form  of  Teflon  until  it  interacts  with  the  arc  itself.  Double  or 
triple  ionization  reduces  the  overall  efficiency  of  the  thruster,  because  of  the  removal  of  more 
electrons  from  one  ion,  instead  of  producing  more  ions.  The  efficiencies  range  from  about  10% 
for  pPPT  up  to  about  90%  for  ion  and  hall  thrusters  [3],  For  electrical  thrusters,  the  energy 
production  capability  of  the  power  plant  limits  the  energy  to  create/accelerate  ions  in  addition  to 
the  associated  electromagnetic  efficiencies  of  acceleration  and  ion  production  itself.  On  a  large 
satellite,  kW  of  electrical  power  may  be  available  in  the  form  of  Radio-Isotope  Thermal 
Generators  (RTG)  or  solar  cells.  However,  for  smaller  satellites  like  CubeSats,  the  power 
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generation  may  only  be  in  the  range  of  watts  or  hundreds  of  watts.  This  low  available  power  is 
one  of  the  main  reasons  pPPT  have  become  an  attractive  option. 

The  Isp  of  PPTs  and  pPPTs  is  typically  around  800-1000s  [3],  Though  they  have  a  lower 
Isp  than  their  electrostatic  or  hall  thruster  counterpart,  their  simple  design  along  with  lower 
power  requirements  make  up  for  any  disadvantage  they  may  have  in  Isp  for  nano  and  pico 
satellites.  Several  spacecraft  launched  commercially,  or  as  part  of  a  tech  demonstration 
employed  PPTs  and  pPPTs  propulsion  devices  onboard. 

History  and  Use  of  PPTs  and  pPPTs  on  small  satellites 

The  devices  are  perfect  for  nano  (1-10  kg)  and  pico-satellites  (0.1-1  kg).  Many  of  these 
satellites  require  thruster  systems  with  high  performance  (Isp)  and  lower  power  consumption. 
Because  these  satellites  are  very  small,  the  perturbations  caused  by  the  impulse  of  chemical 
systems  can  often  be  significant  in  comparison  to  the  satellites  inertia.  The  larger  thrust  to  mass 
ratio  for  smaller  satellites  can  effectively  change  the  orbit  parameters  of  that  satellite  enough  to 
interfere  with  intended  orientation  of  the  satellite  and  would  be  most  notable  in  attitude  changes 
for  the  satellite  where  it  may  possibly  spin  in  its  orbit.  However,  pPPT  are  perfect  for  this  type 
of  satellite  because  the  impulse  bits  are  in  the  range  of  10  to  200  pN-s. 

The  first  PPT  was  flown  on  the  Zond-2  spacecraft  in  1964  launched  by  the  former  Soviet 
Union  [11],  The  first  US  launched  PPT  was  onboard  the  LES-6  mission  in  1968.  In  1976,  the 
LES-8  &  9  satellites  tested  PPTs  for  use  as  station-keeping  devices.  PPTs  have  also  found 
application  for  drag  compensation  for  satellites  typically  below  2000  km  in  altitude,  (above  200 
km  other  influences  dominate,  e.g.  solar  wind,  third  body  effects,  and  photonic  pressure)  [12], 
The  PPT  and  pPPT  were  designed  to  fill  a  niche  for  the  small  satellite  enterprise, 
catering  to  nano  and  pico-satellites.  The  University  of  Washington’s  Dawgstar  satellite  used 
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PPTs  in  formation  flying  experiments.  The  United  States  Air  Force  Academy’s  FalconSat  III 
satellite  used  three  pPPTs  for  two-axis  orbit  stabilization  [13].  As  the  demand  for  smaller  and 
more  precise  propulsion  systems  increase,  pPPT  research  will  continue  to  be  very  valuable  to  the 
commercial/govemment  satellite  ventures. 

II.  Literature  Review 

The  purpose  of  this  chapter  was  to  provide  relevant  operating  principles  and  results  of 
research  to  date  completed  on  the  PPT  and  pPPT  devices  and  implications  to  the  research  topic 
of  this  paper.  The  chapter  begins  by  examining  the  physical  principle  and  physics  behind  the 
operation  of  the  PPT  and  pPPT.  What  is  discussed  next  are  the  relevant  equations  needed  to 
analyze  the  circuitry  of  the  PPT  and  pPPT.  Afterwards,  the  GG-pPPT  circuit  boards  are 
introduced.  From  there,  the  efficiency  of  the  PPT  is  discussed,  and  finally,  the  types  of  failure 
modes  encountered  in  previous  research  is  presented.  Let  us  look  at  the  operation  of  the  PPT 
first  and  then  the  pPPT. 

Pulsed  Plasma  Thruster  and  Micro  Pulsed  Plasma  Thruster 

The  Pulsed  Plasma  Thruster  or  PPT  is  the  predecessor  to  the  pPPT.  It  is  a  larger  device 
and  slightly  more  complicated,  requiring  the  use  of  an  igniter  and  spring  to  maintain  the 
Teflon™  against  the  location  where  the  arc  is  likely  to  form.  The  thermal  ablation  caused  by  the 
arc  does  not  decrease  because  the  Teflon™  has  receded  as  it  is  consumed  to  produce  thrust 
(Figure  4).  The  PPT  being  larger  requires  more  power,  but  is  also  capable  of  producing  thrust  in 
the  pN  to  mN  range.  The  larger  propellant  face  area  exposed  to  the  arc  and  plasma  sheath 
created  in  the  process  of  ionization  results  in  the  higher  thrust  levels.. 

The  pPPT  developed  by  AFRL  at  Edwards  AFB,  is  finding  promising  uses  on  state-of-the  art, 
smaller  satellites.  The  limitation  onboard  energy  production  is  the  biggest  constraint  of 
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propulsion  systems  on  these  spacecraft.  The  pPPT  is  a  simple,  robust  way  to  perform  station 
keeping  maneuvers  and  precise  pointing  movements.  The  pPPT  thruster  shown  in  Figure  5  is  a 
two-electrode  pPPT,  and  was  developed  prior  to  the  three-electrode  pPPT.  The  addition  of  the 
third  electrode  changes  the  operation  from  the  two-electrode  design  as  shown  below  in 
Figure  6.  The  inner  electrode  functions  much  like  the  spark  initiator  for  the  large  scale 
PPT.  This  configuration  has  several  advantages  over  the  two-electrode  design  in  terms  of  break 
down  voltage  and  power  requirements  are  concerned  requiring  less  of  both  to  operate. 


Main  plasma 
sheath 


Figure  6.  Three  electrode  pPPT  operation 

Figure  7,  (a)  shows  the  electrical  connection  for  a  modern  three  electrode  pPPT  with 
electrical  connections  pulsed  power  unit  (PPU)  [14].  The  trigger  capacitor  provides  just  enough 
energy  to  start  seed  plasma  as  the  arc  ablates  and  ionizes  the  inner  annulus  of  Teflon™.  The 
voltage  potential  of  the  trigger  capacitor  must  be  very  high,  in  order  to  overcome  the  dielectric 
breakdown  rating  of  the  Teflon™  surface.  The  electrical  potential  between  the  trigger  capacitor 
and  ground  is  around  nine  kilovolts.  The  impedance  due  to  the  resistor  in  the  trigger  circuit 
controls  the  current  in  the  trigger  circuit,  and  prevents  discharge  of  the  main  from  occurring  too 
early  or  too  late.  This  discharge  then  collapses  the  circuit,  and  the  main  capacitor  discharges, 
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producing  an  arc  that  jumps  from  the  inner  electrode  to  the  outer  electrode  or  anode,  ablating  and 
ionizing  the  outer  ring  of  Teflon™  This  is  the  portion  of  propellant  producing  the  most  thrust. 
The  electrodes  are  made  of  copper,  because  of  its  malleability,  electrical  conductivity  and  heat 
dissipation  properties. 


Outer 

Electrode 


Center  electrope*J 


Main  Capacitor 


(a)  Three-electrode  pPPT,  based  on  Spanjers  [14] 


(b)  Three-electrode  pPPT,  taken  from  Selstrom  [13] 

Figure  7.  3-electrode  pPPT  electrical  connections  and  construction 

A  picture  taken  by  Captain  Jeremy  Selstrom  [13]  shown  in  (b)  of  Figure  7  shows  a  close 

up  of  the  pPPT  and  the  approximate  dimensions  of  the  type  used  in  his  research. 
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Previous  pPPTs  were  originally  of  the  two-electrode  design  as  shown  in  figure  6.  More 
recent  designs  use  a  three-electrode  configuration.  This  design  has  3  significant  advantages  over 
the  two-electrode  design.  First,  the  energy  of  the  main  discharge  has  less  shot-to-shot  variation 
decreasing  the  probability  of  carbon  accumulation  buildup  on  the  propellant  face  [14].  Second, 
the  seed  ionization  produced  from  the  cathode  and  inner  electrode  greatly  reduces  the  voltage 
required  for  the  main  discharge,  which  produces  the  larger  plasma.  Research  indicated  the 
discharge  voltage  could  be  reduced  from  40kV  to  3kV  for  a  14-in  pPPT  if  a  third  or  intermediate 
electrode  was  designed  into  the  thruster  [14].  This  significantly  reduces  the  design  demands  on 
the  PPU.  Third,  the  design  is  more  robust  to  short  term  increases  in  voltages  required  for  initial 
trigger  discharge  and  therefore  arc  formation. 

Electrical  Circuitry  and  Operation 

The  GG-pPPT  firing  sequence  is  very  similar  to  the  operation  of  a  Gatling  gun  machine 
gun,  or  distributor  cap  in  an  automobile.  In  a  Gatling-Gun,  when  the  barrel  rotates  it  collects  a 
round  and  as  the  round  rotates,  it  comes  into  contact  with  the  firing  pin,  which  strikes  the  back  of 
the  round  firing  the  bullet.  As  the  barrel  continues  to  rotate,  the  empty  shell  is  ejected  and 
another  is  fed  into  the  barrel.  This  design  allows  multiple  rounds  to  be  fired  from  the  gun  very 
quickly,  and  still  prevent  the  barrel  from  getting  to  hot.  The  earliest  designs  required  manual 
cranking  to  operate.  Modem  versions  are  powered  with  electric  motors  and  are  gas/liquid  cooled 
because  of  the  extremely  fast  firing  rate. 

A  distributor  cap  in  an  automobile  functions  in  a  similar  way.  A  timing  belt  is  attached 
to  a  pulley  that  rotates  the  shaft  that  is  timed  to  provide  a  spark  to  the  specific  cylinder  1-8  for 
example  in  a  V-8  engine  when  the  piston  is  near  coming  close  to  top  dead  center  (TDC).  As  the 
piston  nears  TDC  the  points  and  condenser,  which  is  just  a  capacitor  provides  a  burst  of  energy 


17 


to  the  specific  spark  plug  in  the  amount  of  about  lOOmJ  of  energy.  This  energy  is  enough  to 
ignite  the  gasoline-air  mixture  vapor  causing  rapid  combustion  and  expansion  forcing  the 
cylinder  back  down  and  providing  power  to  the  wheels  via  a  transmission  and  drive  shaft. 

The  first  step  in  operation  of  a  pPPT  is  to  produce  the  high  voltage  potentials  necessary 
for  arc  formation.  This  voltage  potential  is  accomplished  through  a  resistor-inductor-capacitor 
circuit  or  RLC  circuit.  A  basic  schematic  is  shown  below  in  Figure  8;  the  diode  prevents  the 
current  from  flowing  backwards.  The  circuit  allows  the  use  of  smaller  voltages  in  the  5-28V 
range  for  which  a  small  spacecraft  or  satellite  power  bus  provides.  This  energy  is  stored  in  a 
series  of  capacitors  allowing  a  large  charge  accumulation  to  build  up,  and  therefore  a  large 
voltage  potential  to  be  produced  as  each  switch  in  the  circuit  closes.  When  the  last  switch  closes, 
the  total  accumulated  charge  is  discharged  over  a  period  of  about  20ps  [1],  This  short  discharge 
time  results  in  a  very  large  current  flow  (on  the  order  of  kilo- Amps)  and  a  current  density  of 
MA/m  .  The  current  flow,  J  (Amps)  through  a  capacitor,  resistor,  and  inductor  is  related  to  the 
voltage  via  the  following  equations: 

J=Cd^  (Capacitor)  (1) 

Where  C  is  the  capacitance  in  Farads  dV/dt  is  the  change  in  voltage  per  unit  time. 

J  =  J  ““dt  (Inductor)  (2) 

V(t)  is  the  voltage  at  time  t,  L  is  inductance  in  Henries,  and  t0  and  t  are  the  limits  of  integration 

J  =  (Resistor)  (3) 

R  is  the  resistance  in  ohms.  The  capacitor  and  Inductor  are  energy  storage  elements,  where  the 
resistor  is  a  purely  dissipative  element.  Using  Kirchhoff  s  voltage  laws  and  Thevenin’s  theorem, 
the  voltage  and  current  flow  across  any  element  or  node  in  the  circuit  shown  in  Figure  8.  The 
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differential  equation  that  governs  an  RLC  circuit  can  be  written  in  the  form  of  a  second  order 
homogeneous  differential  equation  [15]: 


L.Q+R.Q+—0 


(4) 


where  Ro  is  the  overall  resistance  of  the  circuit,  with  initial  conditions  Q=Q0,  J0=— Q(0)=0  at 

4  Ln 

t=0.  The  solution  has  two  forms.  For  C<— j-,  there  is  a  damped  oscillatory  behavior,  while  for 

Ro 

4Ln 

C>— ,  the  current  pulse  is  overdamped.  In  general  for  the  under-damped  case,  [15] 

R0 


Q= — -p^=e  2L°  sin(cot+8) 

•  Q 
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For  the  over-damped  case,  [15] 
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The  shorter  the  fall  time  of  the  current,  the  less  propellant  used,  resulting  in  higher 
discharge  densities  allowing  more  ionization  of  the  propellant,  and  producing  more  thrust  and 
higher  Isp..  A  RLC  circuit  with  capacitors  in  parallel  provide  the  momentary,  instantaneous 
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power  required  for  propellant  ablation  and  ionization  as  shown  in  Figure  8.  The  circuit  operates 
in  the  following  way.  Initially,  the  first  switch  is  closed  allowing  the  voltage  to  build  in  the  first 
capacitor.  Then,  the  second  switch  is  closed  allowing  the  voltage  to  build  in  the  second  capacitor 
through  the  inductors,  storing  energy  in  the  form  of  a  magnetic  field,  and  does  not  allow  an 
instantaneous  changes  in  current.  Inductors  are  often  used  in  circuits  as  chokes  to  reduced 
currents  produced  be  by  AC  voltage  sources.  These  components  work  together  to  produce  a  very 
large  voltage  with  the  high  currents  necessary  for  ablation  and  ionization. 


LI  L2  L3 


In  the  GG-pPPT,  a  high-voltage  (HV)  10  kV  multiplier  provided  the  necessary  voltage 
for  firing  each  stick  for  operation.  The  main  capacitor,  high  voltage  switch  board,  external  I/O 
board,  processor  board,  power  input  board,  and  high  voltage  generation  board  are  the  major 
components  of  the  GG-pPPT  module.  The  bundles  of  pPPTs,  arranged  in  a  GG  pattern,  is 
shown  in  (a)  of  Figure  9  and  the  component  part  in  (b)  of  Figure  9.  All  these  components  work 
together  to  allow  sequential  firing  of  each  pPPT  at  various  frequencies. 
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Power  Input  board 

HV  fonorudon 
board 


ExL  VO  board 


Processor  Board 


(a)  Associated  stick  number  and  arrangement  of  (b)  Various  boards  of  GG-pPPT  module  required  for 
pPPT  sticks  Gatling-Gun  Configuration  operation 

Figure  9.  Gatling-Gun  configuration  and  module 

Physical  principles  and  Theory 

The  physical  principle  of  operation  for  all  thrusters  starts  with  momentum  conservation. 
For  electric  thrusters,  the  Lorenz  force  describes  the  acceleration  mechanism  [3]. 


d  1  1  _ y  >  y 

m'"df=q‘<E+u‘XB)+^‘p* 


(7) 


where  m,  is  the  mass  of  the  ion,  du;/dt  is  the  acceleration  of  m,  q,  is  the  particle  charge 
(Coulombs),  E  (V/m)  is  the  applied  electric  field,  B  is  the  magnetic  field  (T),  and  P,k  is  the 
collision  force  per  particle.  However,  we  can  simplify  the  momentum  conservation  for  pPPTs  to 


just  the  contributions  from  the  Lorentz  force.  The  electrostatic  force,  qEdoes  not  provide  the 


majority  of  the  acceleration  of  the  ions  generated  by  the  arc.  The  magnetic  component  qux  B 
imparts  a  significant  force  on  the  plasma,  accelerating  the  ionized  mass  to  speeds  exceeding 
10,000  m/s  for  the  power  and  dimensions  of  these  pPPTs. 

F=q(E+U;XB) 


(8) 


The  effective  exit  velocity  of  the  ionized  Teflon™  is  calculated  in  terms  of  the  radius  of 
the  anode,  cathode,  ablated  propellant  mass  flow  rate,  thrust,  efficiency  and  current  (J)  through 
the  following  equation  [16]: 


ue=— r=— -^-ln 
r|m  r)  471 


R„ 

vRc , 
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m 


(9) 


This  equation  factors  in  the  efficiency  associated  with  the  self-field  magnetic  field 
produced  when  the  Teflon™  interacts  with  the  plasma,  but  is  the  same  equation  as  described  by 
Jahn  with  efficiency  factored  in  [15],  The  term  go,  is  the  permeability  of  free  space  47i  xlO"7 
H/meter  or  N/A2,  J  is  the  current  (Amps),  Ra  is  the  radius  of  the  anode  (m),  Rc  is  the  radius  of  the 
cathode  (m),  m  is  the  mass  flow  rate  (kg/s),  rj  is  the  efficiency,  and  F  is  the  thrust  (N).  Once  the 
exit  velocity  is  known,  Isp  is  determined  through  the  equation: 


Isp  = 


So 


(10) 


, where  g  is  the  acceleration  due  to  gravity  (9.81m/s),  Ue  is  the  exit  velocity  of  the  ions. 

Efficiency 

The  electro-thermal  ablation  process  of  the  propellant  inherently  includes  significant  loss, 
resulting  in  PPTs  only  able  to  achieve  10-20%  efficiency  [3].  The  total  efficiency  written  in 
terms  of  identified  losses  is  [17]: 

ho  =  hppuhtranshsbhhea.hf  (1  1) 


Where  r;0is  the  overall  efficiency,  r]ppu\s  the  power  processing  unit  efficiency,  r)trans is 
the  transfer  efficiency,  ri  Shis  the  sheath  efficiency,  rjheatis  the  heat  loss  efficiency  and  77^  is  the 
frozen  flow  efficiency.  The  r]sh,  r]heat,r]f  efficiencies  are  often  expressed  as  the  acceleration 
efficiency,  r|a  for  the  thruster.  This  term  combined  with  qtrans  is  thruster  efficiency.  The  diagram 
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below  shows  the  relationship  more  clearly.  Total  thruster  efficiency  can  also  be  determined 
based  on  the  exhaust  or  jet  power  by  measuring  thrust  and  the  input  power  supplied.  However, 
multiple  losses  including  waste  heat  generated  by  the  capacitor  and  electronics  are  hidden  in  this 
approach  and  are  only  attributed  to  the  PPU  efficiency.  The  loss  processes  and  associated 
efficiency  statistics  are  graphically  shown  in  Figure  10. 


Figure  10.  Process  efficiencies,  based  on  Bluhm  [17] 

The  exhaust  kinetic  energy  has  two  major  parts,  one  fast  moving  ionized  mass,  m+,  and 


the  other  slower  non-ionized,  neutral  mass,  mn,  for  the  thruster  efficiency  [13]: 
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Where  m+  and  m  are  the  masses  of  the  ions  (F+  and  C+)  and  neutral  particles  in  kilograms 
respectively,  u+  and  u  are  the  exhaust  speeds  of  the  ions  (m/s)  and  neutral  particles,  and  E0  is  the 
input  energy  into  the  system.  Alternatively,  this  can  be  written  in  terms  of  parameters  readily 
measured  from  an  experiment  such  as  thrust,  capacitance  and  mass  flow  rate,  and 
firing/discharge  frequency  [15]. 


F2 

%=- -  (13) 

mvCV2 

where  m  is  the  mass  flow  rate  in  kg/s,  v  is  the  frequency  of  the  discharge  in  Hz,  C  is  the 
capacitance  in  Farads,  and  V  is  the  voltage  supplied  to  the  discharge  capacitor. 

The  total  impulse  of  the  thruster  is  a  combination  of  the  impulse  produced  from  the  ions 
and  neutrals  at  their  respective  velocities. 


:bit=mions+m„eu 

(14) 

I+=m++  m+ 

(15) 

In=m.+  mn 

(16) 

m+ 

— ^=® 

(17) 

The  ratio  of  the  ionized  particles  to  the  neutral  particles  is  designated  by  ®,  important  in 
the  efficiency  of  the  thrusters.  For  different  materials  and  power  levels,  the  ionization  fraction 
varies  considerably  between  zero  and  one.  Multiple  ionization  is  also  possible  for  the  same 
atomic  structure,  where  atoms  can  be  doubly  or  even  triply  ionized  during  their  interaction  with 
the  electric  field.  This  last  term  assumes  a  normal  or  Maxwell  distribution,  which  is  not  a  good 
assumption  for  neutrals  and  ions.  To  do  this,  the  individual  efficiencies  contributing  to  the 
overall  efficiencies  must  be  solved  for.  The  t|ppu  can  be  determined  from  the  energy  transferred 
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to  the  capacitors  from  the  input  power.  This  efficiency  is  normally  between  80-93%  [13].  To 
determine  the  rjtrans,  the  impedance  of  the  PPU  and  the  total  circuit,  Rtotai  must  be  known  and 
from  this  the  efficiency  is  calculated  by  taking  the  ratio  of  the  two  and  subtracting  from  one. 

Ik,,'!1521  (18) 

total 


The  efficiency  of  the  plasma  sheath  can  be  determined  from  the  following  equation, 

Vsh 


hsh=l- 


VD1 


(19) 


Where  Vsh  is  the  sheath  voltage  drop  and  Vppt  is  the  voltage  supplied  to  the  PPT.  The 
heat  efficiency  is  determined  by  looking  at  the  amount  of  heat  generated  from  the  arc  to  the  heat 
lost  that  does  not  result  in  ablation  of  the  propellant, 

'l„..1=l-|t2i  (20) 

E„r„ 


Where  Eioss  represents  the  energy  lost  that  does  not  produce  ablation  and  Earc  is  the 

1 

energy  in  the  arc,  which  can  be  determined  from  the  energy  of  the  capacitors,  —  CV  .  The 


frozen  flow  efficiency,  rjf,  can  be  determined  from  the  speed  of  sound  and  the  total  enthalpy  of 
the  ions,  which  is  determined  from  a  lookup  table  for  the  material.  The  equation  is, 


Tlhe, 


u 


(21) 


Research  performed  by  Selstrom  (USAF)  while  at  the  Air  Force  Institute  of  Technology 
shows  the  effect  of  thruster  efficiency  and  lifetime  as  the  thruster  nears  End  of  Life  (EOL)  [13]. 
As  the  thruster  operated,  the  thrust  and  efficiency  decreased  as  shown  in  Table  1 .  The  reduction 
in  performance  is  due  to  the  larger  carbon  micro-particles  accumulating  between  discharges. 

The  thrust  is  relatively  constant  as  shown  in  Table  2  and  Figure  1 1 .  The  decreased  efficiency  but 
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constant  thrust  would  mean  that  more  propellant  is  coming  off  the  thruster  in  the  form  micro¬ 
particulates  instead  of  the  ionized  atoms  that  are  accelerated  by  the  Lorentz  force. 


Table  1.  Efficiency  as  function  of  firings 


Pulses 

1000 
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10000 

30000 

Efficiency 

(%) 

68.63 

48.19 

48.37 

29.35 

Table  2.  Thrust  as  a  function  of  firings 
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Deviation 

14.20 
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Figure  11.  Thrust  variation  with  firings 


The  pPPT  undergoing  testing  stopped  functioning  at  about  36,000  pulses  due  to  damage 
to  the  working  surface  of  the  thruster.  Since  the  thrust  is  nearly  constant  after  the  initial  phase  of 
operation  and  the  efficiency  decreases,  the  conclusion  of  more  mass  ejected  from  the  thruster  is 
in  the  form  of  neutral  particles  is  supported.  These  particulates  do  provide  a  momentum 
exchange  as  they  leave,  but  are  due  to  their  larger  mass  and  slower  velocity  with  Isp  in  the  range 
of  25-35  s,  instead  of  the  800- 1000s  resulting  from  ionization  and  expulsion  via  the  Lorentz 
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force.  The  premature  failing  of  the  pPPT  stick,  illustrated  one  of  the  major  failure  modes  of  this 
type  of  thruster.  This  mode  is  the  deposit  of  carbon  to  the  propellant  surface  during  thruster 
operation. 

Failure  Modes 

Particulate  Ejection  and  Carbonization 

One  of  the  major  factors  affecting  the  efficiency  is  the  late-time  vaporization  and  carbon 
build  up  on  the  pPPT  as  the  surface  cools  between  pulses  [14].  The  gas  products  from  thermal 
ablation  form  a  carbon  layer  requiring  a  larger  potential  to  pulse.  Contamination  studies  of 
pPPT  by  Kesenek  [18]  reveal  the  carbon  deposits  in  detail  as  shown  in  (a)-(c)  of  Figure  12. 


(a)  20X  magnification  Scanning  Electron  (b)  59X  magnification,  flaky  patches  of 

Microscope  (SEM)  image  of  a  pPPT  stick,  carbon  deposit  are  clearly  visible. 


(c)  1574X  magnification  showing  bumpy  ‘granules’  of  carbon  on  Teflon™  surface 
Figure  12.  Carbon  Buildup,  taken  from  Kesenek  [18] 
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Experiments  and  analysis  conducted  by  Keidar  et  al,  indicates  the  relationship  between 
the  energy  supplied  and  the  amount  of  particulate  that  returns  to  the  surface  of  the  Teflon™  in 
the  form  of  carbon  deposits  [1],  If  the  current  density  is  low,  the  higher  the  ratio  of  returning 
particle  flux  to  ablated  flux  of  Teflon™,  shown  graphically  below  in  Figure  13.  As  evident  from 
the  graph,  the  higher  the  current  density  results  in  more  of  the  ablated  Teflon™  being  ionized 
and  results  in  a  lower  returning  particle  to  ablated  particle  ratio. 


Figure  13.  Particulate/ Ablation  ratio  as  function  of  current  density,  taken  from  Keidar  [1] 

Uneven  Electrode  Erosion 

The  GG-pPPT  can  also  fail  through  the  improper  erosion  of  the  electrodes.  If  the 
electrodes  are  cut  with  a  tool  leaving  ragged  edges  on  the  surfaces  of  the  electrodes,  the  charge 
densities  will  accumulate  on  these  surfaces  because  the  electric  field  will  concentrate  at  these 
micro-protrusions.  The  uneven  erosion  results  in  the  arc  being  to  too  far  from  the  Teflon 
propellant  to  cause  sufficient  ablation  through  radiation  and  convective  processes.  If  the 
electrodes  erode  slower  than  the  propellant,  the  arc  formation  will  be  too  high  above  the  surface 
to  ablate  propellant  effectively,  through  the  electro-thermal  process,  and  will  result  in 
significantly  reduced  thrust  or  will  not  operate.  If  the  electrodes  recede  faster  than  the 


28 


propellant,  the  cathode  may  be  unable  to  see  the  anode  and  the  result  is  poor  operation  and/or 
charring  of  the  propellant,  because  of  the  large  power  required  to  jump  the  gap.  A  diagram  is 
shown  below  in  Figure  14.  Another  failure  mode  is  due  to  uneven  arc  formation  or  ‘spokes’  that 
dig  a  channel  resulting  in  uneven  propellant  utilization  and  ionization. 

Electrodes 


Non-Uniform  Discharge  and  Temperature  Effects 

This  is  another  aspect  of  pPPT  operation  that  is  different  from  the  large  PPT.  The  larger 
PPTs  were  of  rectangular  cross  section,  so  there  was  not  a  point  concentration  of  the  electrical 
current  as  there  is  in  the  coaxial  pPPT.  The  temperature  of  the  propellant  increases  significantly 
where  arc  convergence  is  larger.  These  arc  spokes  are  seen  if  the  energy  of  the  pulse  is  too  small 
[19].  The  small  energy  results  in  non-uniform  arc  distribution,  because  the  current  will  flow 
along  the  path  that  offers  the  least  electrical  resistance,  very  similar  to  the  lightning  strikes  on 
earth  during  severe  thunderstorms.  The  pattern  can  vary  in  both  the  radial  and  azimuthal 
direction  [19].  The  non-uniformity  results  in  higher  propellant  ablation  temperatures  of  the 
propellant  and  affects  efficiency  of  the  thruster. 

The  propellant  utilization  efficiency  suffers  directly  as  a  result  of  this  uneven  propellant 
ablation  and  heating.  Research  conduct  by  Selstrom  [13]  indicates  uneven  ablation  does  indeed 
occur  for  the  pPPT  and  as  a  result,  the  pPPT  firing  rate  would  vary.  A  channel  was  created  in 
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the  Teflon™  as  a  result  of  arc  formation  in  a  concentrated  area.  The  non-ionized  particles 
emitted  from  the  pPPT  due  to  late-time  ablation  does  contribute  to  thrust,  robbing  the  efficiency 
of  the  thruster,  because  these  particles  are  not  accelerated  to  velocities  of  10-20  km/s  by  the 
Lorentz  force.  These  particles  are  more  massive  than  the  ions,  but  are  only  moving  at  sonic 
speeds.  They  do  not  contribute  significantly  to  the  high  Isp  as  shown  in  Table  3. 


Table  3.  Propellant  Velocity  versus  Energy 


Input  voltage  (Energy) 

4  V  (2  J) 

5  V  (2.7  J) 

6  V  (3.4  J) 

7  V  (4.0  J) 

8  V  (  4.4  J) 

Ave  axial  velocity  (m/s) 

215 

256 

345 

270 

258 

Isp  (sec) 

22 

26 

35 

28 

26 

Capt  Seo  Myeongkyo  [20]  in  2007  captured  some  images  of  the  particles  leaving  the 
three  electrode  pPPT  to  characterize  the  exhaust  plume.  The  velocity  of  the  particles  emitted 
ranged  from  100-8076  m/s  with  77%  of  all  particles  collected  traveling  below  300  m/s  as  shown 
in  Figure  15. 


Condition  & 
Velocity 

Image:  X-Z  plane  and  Y-Z  plane 

4V 

8,077  m/s 

30jns 
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1001  m/s 

- - ^ —  J 

Figure  15.  Fast  particles  from  pPPT,  taken  from  Myeongkyo  [20] 

The  low  efficiency  of  the  pPPT  can  be  attributed  to  the  poor  thermal  ablation  process  and 
non-ionization  of  particles  ejected  from  the  device.  Though  it  does  contribute  to  thrust,  the 
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utilization  efficiency  would  be  higher  more  of  the  propellant  is  ionized  before  expulsion  from  the 
thruster.  The  theoretical  efficiency  of  these  devices  could  be  as  high  as  60%  [17].  Current 
performance  is  an  order  of  magnitude  smaller.  The  goal  becomes  how  to  vary  the  input  energy, 
discharge  time  of  the  arc,  and  frequency  of  firing  to  optimize  efficiency  and  propulsion  from  this 
particular  device.  The  GG-pPPT  is  a  collection  of  nine  individual  pPPTs,  and  should  reveal  the 
relationship  between  firing  frequency,  discharge  energy  and  input  power  in  more  detail,  because 
the  electronics  and  their  component  properties  can  be  controlled  via  a  computer  interface  into  the 
unit. 

III.  Methodology 

This  chapter  focuses  on  methodology  and  the  process  of  calibrating  and  setup  of  the 
equipment  in  order  to  conduct  the  research.  We  begin  by  describing  the  facility  and  vacuum 
chamber  where  conduction  of  the  experiment  occurs.  Then,  there  is  discussion  of  the  specific 
equipment  and  procedures  of  connection  and  calibration.  Afterwards,  we  discuss  the  technique 
for  getting  the  center  of  mass  (COM),  moments  of  Inertia  (MOI)  and  physical  parameters  of  the 
torsion  balance.  From  there,  the  research  discusses  how  to  level  the  torsion  balance  and  setup 
the  calibration  electrode  distance.  Then,  we  discuss  the  variables  under  the  operator’s  control. 
Finally,  we  describe  the  planned  test  matrix  and  operating  conditions  under  investigation.  Let  us 
begin  by  describing  the  facility  at  AFIT  where  this  research  took  place. 

The  Geo-orbital  Nano-Thruster  and  Analysis  and  Testing  (GNAT)  Laboratory  at  AFIT 
was  the  primary  testing  facility  for  the  GG-pPPT  research  effort.  The  torsion  balance  required 
calibration  for  testing  of  the  GG-pPPT.  The  available  vacuum  chambers  are  capable  of  reaching 
10~8  torr  [21].  The  chamber  used  in  this  research  was  a  horizontally  mounted  vacuum  chamber 
by  LACO™  Technologies.  The  chamber  pressure  reached  a  low  of  4.7xl0"7  torr,  or  0.47 p torr. 
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The  vacuum  chamber  has  front  and  side  windows  as  well  as  feed-through  (s)  for  high  voltage 
lines  and  BNC  connections  for  instrumentation,  and  test  equipment.  It  is  a  cylinder  roughly  0.61 
m  in  diameter  and  0.61  m  in  depth  with  a  volume  of  0.178  m3.  The  chamber  required  little  input 
from  the  user,  and  operated  using  a  LabView™  interface.  In  order  to  ensure  on-orbit  like 
conditions,  testing  will  occur  below  2x1  O'5  torr,  though  2x1  O'6  torr  is  preferred  according  the 
instruction  manual  for  the  thruster.  Figure  16  shows  a  picture  of  the  chamber  used  with  the  main 
viewport  window  facing  the  camera.  Prior  to  beginning  the  research,  several  key  pieces  of 
equipment  were  required. 

Equipment  and  Connections 

The  most  critical  piece  of  equipment  used  in  this  research  was  the  torsion  balance  system, 
TBS  manufactured  by  Busek™  Space  Propulsion  Company.  Shown  in  Figure  17,  is  the  torsion 
balance.  The  torsion  balance  consisted  of  a  calibration  electrode,  an  active  damping  electrode,  a 
swinging  arm  to  mount  the  thruster  and  its  counterweight,  a  mirror  for  the  laser  displacement 
sensor  or  LDS,  a  mount  for  the  LDS  and  micrometers  to  level  the  balance.  This  device  can 
accurately  measure  pN  forces  produced  by  a  pPPT.  The  LDS  measured  the  linear  travel  of  the 
arm  as  it  rotates,  and  calibrated  according  to  manufacturer  specification  sheets.  The  TBS  also 
has  an  active  damping  system,  but  does  not  require  use  as  long  as  the  permanent  magnet  remains 
in  place.  The  voltage  required  to  operate  the  calibration  electrodes  above  the  noise  floor  present 
in  the  lab  was  set  at  800V  based  on  recommendations  from  Selstrom  in  his  research  [13]. 
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Figure  16.  LACO  Vacuum  chamber  Figure  17.  Torsion  balance 


The  GG-pPPT  shown  in  Figure  18  is  a  nine-stick  array  of  3.175  mm  diameter  3-electrode 
pPPTs.  The  GG-|xPPT  is  a  system  that  requires  a  power  source  capable  of  10-20V  at  four  amps 
or  about  40-80W  peak  in  pulsed  operation.  However,  keep  in  mind  this  power  is  pulsed  and  not 
continuous.  When  the  system  is  idle,  it  only  draws  about  8W.  The  GG-pPPT  comes  with  a  25- 
pin  data  cable  and  a  15-pin  power  cable  that  connects  to  the  device  on  the  top  and  bottom  ports 
respectively  as  shown  in  b  of  Figure  18.  A  re -programming  port  also  exists  on  the  interior  for 
software  changes  shown  in  b  of  Figure  18.  The  re -programming  port  allowed  change  of  the 
discharge  energy  from  the  main  capacitor  in  the  range  from  zero  to  1300V.  However,  to 
reprogram  the  thruster,  removal  of  the  side  panel  is  required,  and  the  high-voltage  wires  had  to 
be  isolated  and  secured  to  prevent  accidental  discharge  of  the  thruster  in  atmospheric  conditions. 
Once  accomplished,  the  thruster  was  ready  to  be  back  into  the  chamber  for  further  testing  at  the 
new  discharge  energy  condition  set  by  the  operator. 
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Figure  18.  External  connections  to  the  GG-|iPPT 
The  next  piece  of  equipment  needed  was  the  stick  selector  module  that  allowed  the  user 
to  change  the  firing  frequency  of  the  thruster  sticks  and  control/enable  which  sticks  are  to  fire. 
The  stick  selector  module  can  run  off  a  nine-volt  battery  for  internal  power,  or  off  external 
power.  Because  of  the  long  wait  times  for  data  processing,  external  power  is  the  preferred  power 
supply  method  for  the  stick  selector  module.  The  center  rotator  knob  controls  the  firing 
frequency  and  the  toggle  switches  enable  or  disable  thruster  sticks.  The  relay  of  status  and  firing 
information  from  the  GG-pPPT  device  to  the  computer  requires  the  RS232  serial  connection.  A 
LED  on  the  side  indicates  power  status.  Shown  in  Figure  19,  is  the  stick-selector  module.  The 
large  bundle  of  wires  needed  for  GG-pPPT  operation  required  attachment  of  a  ‘Cross’  flange  to 
the  top  of  the  vacuum  chamber  in  order  to  feed  all  of  the  connectors  into  the  chamber.  Due  to  its 
length,  the  LDS  tip  is  impossible  to  insert  into  the  chamber  any  way  other  than  straight  down 
from  the  top  of  the  flange  shown  in  Figure  20 
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Figure  19.  Stick  selector  module  Figure  20.  LDS/GG-pPPT  feeds 


An  Agilent  54622D  Mixed  Signal  Oscilloscope  was  used  to  collect  the  raw  data  coming 
from  the  LDS,  and  has  multi-channel  operation  and  digital/analog  signal  acquisition  capability. 
The  voltage  output  of  the  LDS  was  set  in  a  range  where  the  voltage  output  is  a  linear  function  of 
distance  from  the  LDS  tip  to  the  reflecting  mirror.  To  collect  data,  the  output  voltage  of  the  LDS 
required  connection  to  the  oscilloscope  via  a  BNC  connector  and  the  oscilloscope  connected  to  a 
computer.  As  the  thruster  fires,  it  causes  a  small  rotation  in  the  arm  which  shows  up  on  the 
oscilloscope  as  a  change  in  voltage.  Shown  (a)  of  Figure  21  is  the  oscilloscope  along  with  power 
supply(s)  for  the  calibration  electrode  and  the  LDS  box.  The  voltage  changes  due  to  thruster 
firing  are  in  the  mV  range,  so  background  noise  is  a  concern  as  discussed  in  section  four,  results, 
and  analysis.  The  focus  now  shifts  to  collection  of  system  parameters  and  setup. 


(a)  Power  for  Torsion  Balance  and  LDS  (b)  Power  supply  for  GG-pPPT 


Figure  21.  Agilent  54622D  oscilloscope  and  power  supply 
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Center  of  Mass  and  Physical  Properties 
Center  of  Mass 

The  first  step  prior  to  using  the  torsion  balance  was  to  obtain  the  physical  properties  of 
the  thruster.  A  couple  of  ways  to  determine  the  center  of  mass  of  an  object  exist.  The  first  and 
often  simplest  way  is  to  attach  the  object  to  a  slider  able  to  hang  freely  under  the  influence  of 
gravity.  The  hanging  point  is  changed  until  the  center  of  mass  (COM)  is  directly  in  line  with  the 
hanging  wire,  assuming  no  disturbances  are  present,  the  object  should  be  balanced,  with  the  sum 
of  the  moments  and  forces  equal  to  each  other. 

If  the  COM  is  not  at  the  location  of  the  tension  of  hanging  wire,  a  net  moment  will  tend 
to  rotate  the  object  about  its  hanging  point  as  in  Figure  22.  If  the  mass  of  the  object  is  uniformly 
distributed,  then  the  COM  is  usually  the  geometric  center  of  the  object  as  well,  but  does  not  have 
to  be.  If  the  components  of  the  thruster  for  instance  are  more  massive  on  one  side  or  another,  the 
COM  will  shift  from  the  geometric  center.  If  it  is  desired  to  have  the  COM  for  all  three  axes,  the 
thruster  is  hung  from  each  of  its  principle  axis  and  make  adjustments  until  no  moments  in  that 
axis  exist. 
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Figure  22.  COM  determination  by  hanging 
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Another  way  to  measure  the  COM  of  an  object  is  to  place  the  object  on  a  rail  resting  on 
two  knife-edges  fastened  to  a  precision  scale.  First,  the  scales  are  placed  a  measured  distance 
apart.  The  center  of  the  rail  was  placed  at  the  midpoint  distance  between  the  scales.  The  COM 
of  the  rail  was  determined  by  adjusting  the  location  of  the  rail  until  the  scales  read  the  same 
value.  The  thruster  was  placed  on  the  rail  and  moved  back  and  forth  until  both  scales  read  the 
same  value  to  within  a  milligram  of  each  other.  The  change  in  distance  from  the  geometric 


center  of  the  thruster  along  with  the  scale  measurements  allowed  precise  determination  of  the 
COM.  Because  of  the  transformers  and  various  components,  the  COM  was  off  from  the 
geometric  center  by  0.5  mm.  The  knife-edge  setup  used  in  this  research  and  the  precision  digital 
caliber  to  measure  the  dimensions  of  the  thruster  is  shown  in  Figure  23. 
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(a)  Digital  scales  with  knife  edge  and  rail  (b)  Digital  caliper  for  measuring  distances 

Figure  23.  Knife  edge/scales  to  determine  center  of  mass 

The  scales  were  capable  of  reading  to  0.0001  kg,  (0.1  g)  with  a  capacity  of  6  kg.  This 
process  was  repeated  for  each  principle  axis  of  the  thruster.  With  the  COM  of  the  thruster 
known,  the  next  task  was  to  make  physical  measurements  of  the  torsion  balance  itself  and 
determine  the  moment  arms  of  the  thruster  and  counterweight. 

Moment  Arms 

The  next  step  prior  to  using  the  torsion  balance  is  to  characterize  its  physical  dimensions, 
because  these  dimension  were  important  in  determine  the  moment  arms  of  the  thruster, 
counterweight,  and  location  of  the  LDS  sensor  tip  relative  to  the  axis  of  rotation  for  the  torsion 
balance,  which  is  was  designated  zero  for  reference  purposes.  The  length  of  the  arm  of  the 
torsion  balance  was  measured  using  a  combination  of  measuring  tapes  and  the  calipers  shown  in 
Figure  23.  The  length  of  the  torsion  balance  arm  measured  0.476  m  from  end  to  end,  with  the 
center  bar  having  a  thickness  of  25.4  mm.  The  base  plates  of  the  balance,  which  hold  the 
counterweight  and  thruster  are  squares  that  are  38.1  mm.  on  each  side.  The  distance  from  the 
center  of  one  plate  to  the  other  is  therefore  0.400  m. 

Because  of  the  wiring  required  to  connect  the  thruster  to  its  power  supply,  the  radius  of 
the  thruster  is  not  exactly  the  same  as  the  radius  of  rotation  for  the  counter  weight.  They  differ 
by  about  3.17  mm.  The  location  of  the  LDS  sensor  was  found  by  measuring  the  distance  from 
the  end  of  the  balance  arm  to  the  center  of  the  mirror,  accounting  for  the  diameter  of  the  LDS 
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sensor  tip.  This  The  result  is  0.225 12  m  for  the  LDS  sensor  tip  and  a  radius  of  rotation  of  0.200 
m.  for  each  arm  if  the  counterweight  is  located  at  the  center  of  each  mounting  plate.  With  these 
values  known  for  the  balance,  the  next  task  was  to  determine  the  masses  and  moments  of  inertia 
(MOI)  for  both  the  counterweight  and  thruster.  The  moment  of  inertia  of  the  balance  arm  was 
provided  by  Busek™. 

Calculation  of  Mass  and  Moment  of  Inertia  (MOI) 

The  mass  of  the  counterweight  and  the  thruster  were  measured  directly  using  the  same 
precision  scales  used  to  determine  the  center  of  mass  of  the  thruster.  The  next  step  afterwards, 
was  to  figure  out  the  MOI  of  both.  The  mass  of  the  thruster  with  its  stand  was  1.9208  kg;  the 
mass  of  the  counterweight  with  its  base  was  1.921  kg.  These  masses  differ  by  only  0.2  mg.  The 
error  associated  with  the  difference  in  mass  is  only  0.01%  of  the  total  mass  of  each.  The  MOI 
for  both  the  counterweight  and  the  thruster,  can  be  calculated  using  the  parallel  axis  theorem  and 
MOI  formulas  about  the  body  axis  of  the  counterweight  and  thruster.  The  thruster  and  counter 
weight  are  both  rectangular  solids  in  by  design,  so  this  makes  the  calculation  of  the  MOI 
relatively  simple.  The  MOI  of  a  rectangular  solid  about  its  own  axis  of  rotation,  say  the  z-axis, 
which  we  will  call  the  height  of  the  thruster  and  counterweight  is: 

t=^m(w2+d2)  (22) 

The  terms  w  and  d  are  the  width  and  depth  of  the  rectangular  solid  respectively.  Once  Iz 
was  determined,  the  parallel  axis  theorem  was  used  and  the  total  MOI  was  determined  from  the 
equation: 

(23) 

The  term  r  is  the  radius  of  rotation  from  its  body  rotation  axis  to  the  desired  axis  of 
rotation,  and  is  equal  to  the  radius  between  the  center  of  the  torsion  balance  and  the  COM  of  the 
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thruster  or  the  counterweight  in  the  X-Y  plane.  The  moment  MOI  for  the  thruster  using  formulas 
above  was  0.0791  kg-m2  and  0.08037  kg-m2  for  the  counterweight  respectively.  The  MOI  for 

the  balance  arm  is  0.01692  kg-m  provided  by  Busek.  The  total  MOI  is  just  the  sums  of  these 

2 

which  are  equal  to  0.17643  kg-m  . 

Calibration  Pulse  and  Force  Calculation 

The  torsion  balance  does  not  measure  thrust  directly,  but  calculates  it  through  a 
comparison  to  a  calibration  pulse  of  known  voltage  and  a  force  conversion  factor  that  is  a 
property  of  the  geometric  properties  and  spacing  between  the  electrodes.  Assuming  the  spacing 
between  the  electrodes  is  1.0  mm,  the  force  can  be  calculated  using  electromagnetic  principles 
and  modified  equations  for  the  attraction  between  two  parallel,  electrically  conducting  plates.  In 
collecting  data,  a  calibration  pulse  should  be  taken  as  frequently  as  possible  to  obtain  more 
accurate  results  and  calculations.  These  values  can  be  statistically  analyzed  to  see  if  the 
variance,  mode,  or  standard  deviation  changes  significantly  between  one  run  and  the  next.  Once 
the  calibration  pulse  is  completed,  the  calibration  voltage  is  turned  off.  After  some  settling,  the 
thruster  is  started  and  data  collection  can  take  place.  A  calibration  pulse  should  have  an 
exponential  sinusoidal  decay  based  on  the  physical  dynamics  of  the  torsion  balance.  It  is 
essentially  an  oscillator  with  viscous  forces  resisting  the  movement  of  the  balance  arm.  A  good 
calibration  pulse  can  be  seen  on  the  oscilloscope  for  qualitative  assessment  as  in  Figure  24. 
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(a)  Oscilloscope  display  of  calibration  pulse  and 

thruster  response  after  settling.  (b)  MATLAB™  piot  of  calibration  pulse 


Figure  24.  Typical  calibration  pulse 

The  data  collected  for  various  calibration  pulses  and  the  thruster  is  presented  in  the 
analysis  and  results  chapter  of  this  thesis. 

Calibration  of  Laser  Displacement  Sensor  (LDS) 

The  TBS  requires  the  LDS  operation  in  a  region  where  the  voltage  is  a  linear  function  of 
displacement.  The  LDS  is  manufactured  by  Philtec™  and  the  manufacturer  provided  sensitivity 
measurements  and  calibration  curves.  The  near  side  sensitivity  of  the  LDS  is  0.916  mV/pin. 

The  far  side  sensitivity  slope  is  -15.30mV/mil,  and  used  to  relate  the  voltage  change  to  a  change 
in  linear  displacement  of  the  arm  and  mirror  relative  to  the  LDS  sensor.  For  the  sensor,  a  voltage 
range  of  three  to  three  and  a  half  is  used  as  the  operating  region  of  the  LDS,  and  is  where  the 
output  voltage  varies  directly  with  distance  between  reflecting  mirror  and  the  LDS  sensor  tip. 

The  LDS  has  coarse  and  fine  correction  screws  and  is  adjusted  so  the  maximum  voltage 
output  from  the  LDS  is  five  volts.  The  LDS  operates  on  either  the  ‘near’  and  ‘far’  side,  with 
each  calibration  curve  having  different  constants,  so  the  operator  needs  to  know  which  side  of 
the  curve  applies  to  the  system  under  operation.  Once  the  LDS  unit  is  calibrated,  the  next  step  is 
to  determine  the  proper  distance  the  calibration  electrodes  need  to  be  apart. 
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Calibration  of  Electrode  Distance 


The  calibration  electrode  distance  must  be  exactly  1mm  apart  to  use  the  voltage  force 
conversion  formula  in  the  TBS  operations  manual.  If  the  electrodes  move  too  close  together,  and 
make  contact  each  other,  they  will  short  out.  If  shorted  out,  the  electrodes  will  remain  together 
and  calibration  data  cannot  be  taken.  This  occurred  twice  during  the  research,  and  each  time,  the 
chamber  brought  back  to  atmospheric  conditions,  recalibrated,  and  pumped  back  down  to 
vacuum  conditions.  Once  the  electrodes  are  the  correct  distance  apart,  the  force  of  attraction 
between  the  electrodes  is  calculated  using  equation  (24) 

This  equation  is  a  modified  form  of  the  equation  for  two  parallel  conducting  plates  but 
adjusted  due  to  the  different  sizes  of  the  electrodes  and  electric  field  variations  near  the  edge  of 
the  electrodes.  This  set  distance  is  accomplished  in  the  following  way.  First,  the  operator 
verifies  the  LDS  system  is  on  and  calibrated  as  stated  above.  Second,  have  the  spacing  initially 
about  1.0  mm  apart,  and  move  the  calibration  electrode  towards  the  electrode  on  the  TBS  arm 
until  the  LDS  voltage  changes  up  or  down  as  shown  in  Figure  25. 


Figure  25.  Electrode  calibration 
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Once  the  electrodes  make  contact,  there  should  be  an  immediate  vertical  shift  in  the 
oscilloscope  trace.  Back  the  electrode  off  until  this  shift  is  as  small  as  possible.  Third,  using  the 
position  micrometer  turn  the  micrometer  two  full  turns  counterclockwise.  Each  full  turn 
corresponds  to  500  pm,  so  two  full  turns  is  1000  pm  or  1.0  mm.  The  electrodes  are  now  exactly 
1.0  millimeter  apart  in  spacing.  The  formula  to  calculate  the  calibration  force  between  the  set  of 
electrodes  is: 


(1.262xl0"9)(l+138.1*d)V 


cal 


(24) 


The  quantity  d  is  the  distance  between  the  electrodes  (1.0  mm)  if  set  properly  using  the 
procedure  discussed  in  the  operations  manual  above.  Vcai  is  the  applied  calibration  voltage, 
which  was  800V  in  this  case  to  provide  enough  force  to  move  the  balance  arm  above  the  noise 
floor.  This  force  is  inversely  proportional  to  the  square  of  the  distance  between  the  electrodes 
and  proportional  to  the  square  of  the  applied  voltage.  The  non-linearity  of  this  equation  poses 
challenges  if  the  electrodes  move  significantly,  especially  if  it  results  in  moving  outside  the 
linear  response  range  of  the  LDS. 

The  electrode  distance  must  be  set  prior  to  obtaining  vacuum  conditions,  the  LDS 
response  may  drift  due  to  variations  in  thermal  conditions  and  metallic  expansion/contracting 
inducing  strain.  To  correlate  the  effect  of  thruster  with  the  oscilloscope,  the  sign  convention  and 
direction  must  be  established.  The  sign  convention  of  the  system  is  shown  in  Figure  26.  From 
Figure  26,  if  the  electrodes  are  moving  move  further  away  from  each  other,  then  the  mirror  is 
moving  closer  to  the  LDS  tip  resulting  in  an  increase  in  the  voltage  output  from  the  LDS.  This 
result  indicated  the  LDS  is  operating  on  the  “far”  side  of  its  sensitivity  curve.  If  it  were  on  the 
near  side,  as  the  gap  became  smaller,  the  voltage  output  would  decrease.  The  near  side 
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sensitivity  slope  of  0.916  mV/pin  would  be  ideal  for  sub  micro-Newton  thrusters,  but  proved  too 
difficult  to  use  in  practice. 
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Figure  26.  Sign  convention/leveling  devices-using  far  side  of  LDS 
The  extreme  sensitivity  on  the  near  side  of  the  sensor,  made  that  side  impractical  to  use. 
Use  of  the  near  side  is  desired  for  sensitivity  reasons,  but  was  too  sensitive  and  would  quickly 
move  out  of  the  linear  operating  range  of  the  LDS,  or  would  saturate  and  move  to  the  topside  of 
the  curve,  making  it  useless  for  data  gathering  and  analysis.  The  near  side  of  the  LDS  only  has  a 
useful  range  of  1.9  (0.05  mm)  to  3.7  mil  (0.094mm).  Thermal  drift  and  external  noise  are 
enough  to  move  the  TBS  out  of  this  range.  However,  the  far  side  has  a  useful  range  from  37  mil 
(0.942  mm)  to  149.2  mil  (3.79  mm),  which  should  be  controllable. 

Balancing  and  Leveling  the  Torsion  Balance 

The  LACO™  vacuum  chamber  only  has  a  volume  of  0. 178  m  ,  and  made  it  very 
confined  space  to  work  inside  with  the  torsion  balance,  thruster,  and  counterweight.  A  special 
stand  constructed  out  of  80/20™  aluminum  was  used  to  elevate  the  torsion  balance,  which  was 
mounted  to  an  optical  bench  plate.  The  adjustment  screw  on  the  left  side  of  cross  arm  in  Figure 
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26  was  inaccessible  so  two  more  screws  were  added  to  the  optical  bench  to  allow  leveling  in  2- 
axes.  The  additions  of  these  screws  provided  an  additional  benefit  of  providing  greater  control 
of  the  leveling  procedure  without  having  to  bump  the  counterweight  or  disturb  the  wiring  inside 
the  chamber. 

The  leveling  of  the  balance  should  be  done  with  the  thruster  and  counterweight  installed, 
because  the  weight  of  these  components  will  affect  the  leveling  and  strain  on  the  stand  and  need 
the  level  point  may  change  if  these  mass  effects  are  not  taken  into  account  in  the  baseline 
measurements.  The  counterweight  was  milled  out  of  T6061  aluminum  and  had  dimensions  of 
(127  x  76.2  x76.2)  mm,  cross  sectional  area  of  5.80644  *10'3  m2  and  a  volume  of  7.374*10'4  m3. 
Using  the  volume  and  mass,  the  density  of  the  block  is  calculated  to  be  2604  kg/m  .  The  block 
of  aluminum  also  had  a  7.9375  mm  hole  drilled  into  the  bottom,  so  it  could  be  mounted  on  a 
sheet  metal  base  and  attached  to  the  thruster  mounting  plate.  It  was  recommended  later  the 
counterweight  not  be  fixed  to  compensate  for  any  center  of  gravity  (CG)  error  due  to  moving  the 
thruster  in  and  out  of  the  chamber. 

One  of  the  other  pieces  of  equipment  needed  to  level  the  balance  properly  was  a  pivoting 
mirror.  The  tight  space  of  the  chamber  prohibited  anyone  from  sticking  their  head  in  to  directly 
view  bubble  level  provided  with  the  torsion  balance.  So,  this  mirror  is  required  in  combination 
with  a  laser  and  LED  light  source  to  see  when  the  bubble  was  in  the  center  of  the  black  circle 
indicating  the  system  was  level.  A  device  was  constructed  to  allow  viewing  of  the  bubble  level 
and  the  physical  operation  of  the  thrust  using  two  mirrors  mounted  to  an  80/20™  stand  at  45 
degrees,  effectively  acting  as  a  beam  redirector.  The  stand,  mirror,  and  the  bubble  level  are 
shown  in  a-d  of  Figure  27.  The  picture  in  (a)  of  Figure  27  redirected  the  image  of  the  thruster 
out  of  one  of  the  port  windows  on  the  side  of  the  chamber. 
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The  size  of  the  chamber  severely  limited  pointing  options  for  the  thruster.  If  the  thruster 
is  placed  in  the  front,  it  would  not  be  possible  to  adjust  the  leveling  micrometers  from  the  back 
of  the  vacuum  chamber  without  disturbing  the  test  setup  significantly.  The  length  of  the  balance 
optical  bench  and  the  arm  would  also  not  permit  sideways  mounting  of  the  torsion  balance, 
where  the  thruster  would  face  the  main  viewing  window  on  the  front  of  the  chamber.  The  other 
difficulty  in  repositioning  the  torsion  balance  is  to  avoid  the  wires  going  to  the  GG-pPPT  unit 
from  twisting  too  much.  Any  excessive  twisting  in  the  wires  will  force  them  to  act  like  a 
torsional  spring,  adding  undesired  responses  and  measurements  of  the  dynamics  of  the  system,  if 
severe  enough.  There  is  also  concern  of  damaging  the  fiber  optic  fibers  of  the  LDS  feeding  into 
the  chamber.  Excessive  twisting  could  potentially  result  in  fracture  of  these  fibers  and  degraded 
LDS  performance  and  measurements. 


(a)  Thruster  observation 
stand 


(b)  Pivoting  mirror  to  see  position  of  bubble 
leveler 


(d)  LED  illumination 

Figure  27.  Viewing  mirrors  and  bubble  leveler 
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The  image  in  (b)  of  Figure  27  allowed  the  leveler  to  be  viewed  from  the  top,  because  the 
mirror  placed  at  a  45-60  degree  angle  to  the  vertical  provided  a  top  view  of  the  bubble  leveler, 
which  eliminating  the  need  for  an  individual  to  stick  their  head  into  the  chamber  and  risk 
perturbing  the  system.  With  the  balance  level,  connections  made,  and  calibration  of  LDS  and 
torsion  balance,  we  now  discuss  the  input  parameters  that  the  operator  can  vary,  and  the  test 
point  matrix. 

Independent  Control  Variables 

The  four  parameters  controlled  in  the  operator  of  the  GG-pPPT  are  discharge  energy, 
firing  order,  pulse  frequency,  and  input  voltage.  The  discharge  energy  of  the  arc  is  crucial  due  to 
operation  and  potential  carbon  deposition.  The  firing  order  can  reveal  the  impact  of  buildup  and 
cross  contamination  as  well  as  potential  asymmetric  thrust  from  one  stick  to  another.  The  input 
frequency  can  vary  from  0.5  to  2.0  Hz.  The  higher  frequency  provides  the  ability  to  achieve  a 
pseudo-constant  thrust  level.  The  variation  in  input  voltage  allows  determination  of  the 
minimum  power  required  for  thruster  operation. 

The  GG-pPPT  system  is  software  controlled  with  an  RS232  serial  connection.  The 
operator  of  the  thruster  can  see  the  status  of  the  individual  pPPTs  on  the  Cubloc™  software 
provided  by  Busek™.  Once  all  the  connections  are  made,  the  operator  provides  power  to  the 
GG-pPPT  unit  and  opens  up  the  debugger  tool  in  the  Cubloc™  software  that  came  with  the  unit. 
The  operator  should  see  a  Gatling-gun  pattern  of  ‘X’s,  indicating  those  sticks  are  not  enabled. 
When  enabled,  the  ‘X’s  will  turn  into  a  number,  indicating  that  stick  is  active  and  ready  to  fire. 
The  operator  then  determines  the  sequential  firing  order  of  the  sticks  by  turning  on  the 
appropriate  toggle  switches  numbered  one  thru  nine  on  the  stick  selector  module. 
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The  operator  then  turns  the  frequency  knob  to  a  value  between  a  one-half  hertz  and  two 
hertz.  The  sticks  are  now  firing  and  should  show  up  as  voltage  changes  on  the  oscilloscope  as  in 
(a)  of  Figure  24.  A  complete  setup  of  the  system  is  shown  in  Figure  28.  The  necessary 
calibration  voltage  needed,  is  provided  using  the  EMCO™  DC-DC  power  converter  and  a 
Hewlett-Packard™  200W  power  supply.  The  800V  voltage  was  routed  into  the  chamber 
separate  from  the  other  connections  to  prevent  the  possibility  of  arcing  or  short-circuiting  to  a 
component  in  the  thruster.  The  power  to  electrodes  was  supplied  only  when  a  calibration  pulse 
needed  to  be  taken,  for  instance  when  the  thruster  was  reprogrammed  and  placed  back  into  the 
vacuum  chamber  for  further  testing. 


Data  TX/RX 


48 


To  determine  the  effect  of  frequency,  discharge  energy,  and  firing  pattern  on  the 
performance  of  the  thruster,  a  test  matrix  was  required.  The  test  conditions  are  repeated  enough 
times  so  bias  and  experimental  error  can  be  determined.  This  thruster  configuration  and  module 
being  new,  very  little  research  into  its  operation  and  characterization  exist. 

The  testing  to  date  on  this  thruster  was  simply  to  determine  the  impulse  bit  or  Ibit  for  the 
thruster,  and  the  average  mass  loss  per  discharge.  The  matrix  of  planned  test  conditions  is  shown 
in  Table  4.  This  combination  of  parameters  provided  192  possible  test  conditions,  more  than 
sufficient  to  provide  the  needed  data  sample  size  for  statistical  data  analysis  and  comparison. 

Table  4.  Test  Conditions  for  various  input  voltage,  frequency  and  discharge  energy 


Discharge 

energy 

(J) 

E= 

Output 
Voltage  (V) 

=l/2CV2,C=2pF 

Frequency(Hz) 

Stick  pattern 

Voltage  (V) 

1.7 

1303.84 

0.5, 1,1.5, 2 

1,2, 3, 4, 5, 6, 7, 8, 9 

12.5,15,17.5,20 

1.6 

1264.91 

0.5, 1,1.5, 2 

1,2, 3, 4, 5, 6, 7, 8, 9 

12.5,15,17.5,20 

1.5 

1224.74 

0.5, 1,1.5, 2 

1,2, 3, 4, 5, 6, 7, 8, 9 

12.5,15,17.5,20 

1.4 

1183.22 

0.5, 1,1.5, 2 

1,2, 3, 4, 5, 6, 7, 8, 9 

12.5,15,17.5,20 

1.3 

1140.18 

0.5, 1,1.5, 2 

1,2, 3, 4, 5, 6, 7, 8, 9 

12.5,15,17.5,20 

1.2 

1095.45 

0.5, 1,1.5, 2 

1,2, 3, 4, 5, 6, 7, 8, 9 

12.5,15,17.5,20 

Discharge 

Output 

Frequency  (Hz) 

Stick  pattern 

Voltage  (V) 

energy 

Voltage  (V) 

(j> 

1.7 

1303.84 

0.5, 1,1.5, 2 

3,8 

12.5,15,17.5,20 

1.6 

1264.91 

0.5, 1,1.5, 2 

3,8 

12.5,15,17.5,20 

1.5 

1224.74 

0.5, 1,1.5, 2 

3,8 

12.5,15,17.5,20 

1.4 

1183.22 

0.5, 1,1.5, 2 

3,8 

12.5,15,17.5,20 

1.3 

1140.18 

0.5, 1,1.5, 2 

3,8 

12.5,15,17.5,20 

1.2 

1095.45 

0.5, 1,1.5, 2 

3,8 

12.5,15,17.5,20 

After  the  equipment  was  set-up  and  the  necessary  connections  are  made  as  shown  in 
Figure  28,  the  physical  properties  of  the  thruster,  counterweight  and  torsion  balance  was 
measured  and  collected.  Once  complete,  the  equipment  was  calibrated  and  baseline  data  taken  to 
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establish  a  reference  frame  from  which  the  data  collected  was  analyzed.  We  are  now  to  ready  to 
start  collecting  data  from  the  thruster  and  start  by  solving  the  equations  of  motion  (EOM)  for  the 
torsion  balance  system  and  its  thruster  and  counterweight. 

Solution  approach 

The  desired  method  of  modeling  for  the  torsion  balance  system  and  the  thruster  was  to 
model  the  system  as  a  harmonic  oscillator  and  solve  the  equations  of  motion  directly  to  achieve 
predicted  deflection  angles,  theta  in  radians.  Once  this  deflection  angle  is  known,  the  force  and 
hence  thrust  from  the  thruster  can  be  found  by  applying  rotational  dynamics  and  statics  to  the 
system  and  find  the  force  that  produces  the  observed  deflection  with  the  determined  spring 
constants  for  the  entire  system.  This  method  was  iterated  until  the  root  mean  square  error 
(RMSE)  was  minimized  or  matching  was  optimal,  based  on  the  pulse  train  of  the  input  signal  for 
various  thrust  levels  from  0  to  200  pN.  The  formula  is  fairly  simple  and  has  the  form: 


Yn  (X  data-X  model)2 
RMSE  =  J^1=1  ~ - = - — 


(25) 


RMSE  has  the  same  units  as  the  dependent  variables  and  is  the  error  doubles,  the  RMSE 
doubles.  This  direct  relationship  is  not  so  if  using  MSE,  but  it  is  still  widely  used.  The  smaller, 
the  RMSE,  the  better  fit  the  model  is  to  the  data  in  question.  The  solution  approach  can  also  be 
depicted  graphically  and  shown  Figure  29. 
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]  Determine  £  and  con  from  calibration 
via  T  and  tod  and  iterate  until  convergence 

^J=lace  values  of  ^  and  ton  into  State-Space 
representation  of  system 


6 


input  force  values  with  saw-tooth  or  square  waveform 
at  desired  frequencies 


Q  Run  model  and  calculate  RMSE 
between  model  and  data, 
iterate  until  min  error  is  achieved 


Verify  in 
Simulink! 


Figure  29.  Graphical  Solution  Approach 


IV.  Analysis  and  Results 

This  chapter  begins  by  looking  at  the  mathematics  and  physical  properties  of  the  torsion 
balance  system,  and  from  there,  develops  the  equations  of  motion  (EOM)  for  the  torsion  balance 
and  its  relation  to  the  force  provided  by  the  thruster.  Next,  some  plots  for  the  range  of  discharge 
energies  from  1.7  to  1.2  J  comparing  the  data  gathered  against  the  model.  Then,  discussion  of 
power  needed  for  operation  and  the  accumulation  of  carbon  on  the  propellant  face.  Afterwards, 
the  discussion  turns  to  the  variation  of  thrust  with  discharge  energy,  frequency  and  input  power. 
Next,  we  discuss  the  results  of  determining  torque  from  the  firing  of  sticks  three  and  eight  of  the 
GG-pPPT  module.  Finally,  we  discuss  the  various  amounts  and  forms  of  error  present  in  the 
experiment 
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Equations  of  Motion 

The  first  step  in  understanding  the  behavior  of  the  thruster  was  to  solve  the  equations  of 
motion  for  the  torsion  balance  and  to  use  this  solution  to  model  the  behavior  of  the  balance  due 
to  the  thruster  for  various  firing  frequencies  and  discharge  energies.  The  equation  of  motion  of 
the  torsion  balance  was  a  second-order  ordinary  differential  equation,  ODE.  The  solution  of  this 
ODE  required  use  of  state-space  methods  or  Laplace  Transforms  from  linear  systems  theory. 

From  linear  systems  theory,  the  force  is  not  always  a  linear  function  of  the  spring 
constant  if  the  angle  is  very  large,  because  all  springs  have  a  range  where  their  force  is  no  long 
linearly  proportional  to  the  amount  of  deflection.  This  phenomenon  is  because  of  elastic  limits 
of  real  materials.  Fortunately,  the  angles  resulting  from  the  thruster  are  approximately  10'5 
radians,  so  small  angle  approximations  are  valid  for  this  analysis  and  provided  excellent  results. 
The  angle  of  displacement  can  be  determined  from  some  trigonometry  and  data  about  the 
moment  arm  and  LDS  measurements  as  in  Figure  30.  For  angles  on  order  of  tens  to  micro¬ 
radians,  the  inverse  tangent  and  inverse  sine  functions  provided  the  same  result  essentially, 
because  of  the  small  angles  involved. 

rLDS  _tip*sin(0) 

=Y(displacement) 


Figure  30.  Trigonometry  to  determine  angular  displacement 
Determination  of  the  angle  was  accomplished  by  taking  the  inverse  sine  of  the  ratio  of  the 
displacement,  to  the  radius  of  rotation  of  the  LDS  sensor  tip.  The  displacement  was  determined 
by  converting  the  LDS  output  voltage  into  an  equivalent  displacement.  For  example  if  the 
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sensitivity  was  15.30  mV/mil  (602  mV/mm)  and  we  saw  a  1.0  mV  change  from  the  LDS,  then 
the  displacement  is  1/602  mm  or  -1.661 1*10"3  mm  or  1.661 1  pm.  This  relationship  is  shown  in 
Eq.(26) 


The  angle  and  hence  the  displacement  from  the  neutral  position  was  found  through  the 
above  method.  This  assumption  allows  the  equations  of  motion  to  be  modeled  as  a  linear 
system.  With  this  being  the  case,  the  angular  displacement  was  directly  proportional  to  the 
spring  constant  of  the  system.  The  torsion  balance  has  an  internal  spring  constant  of  1 1.6 
mN/radian,  but  the  spring  constant  of  the  wiring  and  connecting  structure  must  be  taken  into 
account  to  get  an  “equivalent”  spring  constant  for  the  system.  This  effect  should  not  affect  the 
overall  thrust  calculation,  but  rather  will  result  in  a  stiffer  spring  constant  and  therefore  smaller 
angular  displacement. 

The  system  model  for  the  balance  and  thruster  is  essentially  a  harmonic  oscillator 
described  by  the  second  order  ODE.  The  parameters  that  affect  the  behavior  are  the  damping 
ratio,  and  the  natural  frequency  of  the  system,  (on.  The  equation  for  the  system  is  [22]: 


e^co  e+co  2e=h!& 

^  n  n  t 

total 


The  first  and  second  derivatives  of  angular  position  are  angular  velocity  and  angular  acceleration 
respectively.  F(t)  is  the  forcing  function  (thruster)  and  Rthmster  is  the  moment  arm  of  rotation,  and 
Total  is  the  total  moment  of  inertia  for  the  torsion  balance  including  the  arm,  thruster  and 
counterweight.  If  the  damping  ratio  is  less  than  one,  the  system  will  have  at  least  one  complete 
oscillation  and  is  under-damped.  If  the  damping  ratio  is  larger  than  one,  the  system  will  have  no 
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oscillations  and  critically  damped.  The  system  will  respond  like  a  first  order  differential  system 
where  it  will  rise  steeply,  but  then  asymptotically  approach  its  equilibrium  or  final  value  at  t-^oo. 
This  system  is  the  angular  analog  to  the  spring  mass  damper  system. 

The  damping  ratio  can  be  determined  by  looking  at  the  calibration  pulse.  The  natural 
frequency,  con  can  be  determined  from  cod,  which  is  the  damped  natural  frequency  of  the  system, 
and  the  only  one  observed  experimentally.  They  are  related  to  each  other  through  the  damping 
ratio,  £,  so  once  one  is  known,  the  other  can  be  found.  The  ability  to  determine  the  quantities 
accurately  along  with  measurement  errors,  system  noise,  and  tolerances  will  dictate  how  closely 
the  model  matches  the  data.  Ziemer  [22]  solved  the  equation  for  an  impulse  input.  The  initial 
conditions  for  model  at  time  equal  to  zero  was  determined  from  the  data  when  the  pulse  first  was 
detected  by  the  oscilloscope  and  Cublock™  software.  The  impulse  solution  to  the  above 
equation  has  the  form: 


0(t)= 


I  *  R  -£* co 

bit  thruster  fl 


I  * 

total 


CO 


d 


*t 

sin((Btj*t+(I>) 


(28) 

The  term  <D  is  the  phase  shift  of  the  signal  in  radians  from  zero.  For  the  system  at  rest, 
the  phase  angle,  ®  should  be  zero  for  the  simplest  case.  The  impulse  bit  or  Ibit  is  determined 
from  the  initial  velocity  at  t=0,  Itotai,  and  radius  of  rotation,  Rthmster  through: 


I 


bit 


0*1 


total 


R 


thruster 


(29) 


The  angular  velocity  at  time  zero,  can  be  evaluated  by  determine  the  change  in  slope  of  the 
angular  displacement  very  close  to  the  time  when  the  input  supplied  force,  F  (t).  Some  error 
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results,  but  if  the  time  step  is  very  small,  it  should  provide  a  good  approximation  for  angular 
velocity  at  t=0.0.  Once  these  quantities  were  determined,  MATLAB™  can  model  the  system 
and  compare  that  result  to  the  data  for  comparison  and  amount  of  agreement  between  theory  and 
experiment.  Laplace  Transforms  are  another  method  for  determining  the  response  of  the  system 
if  the  forcing  function  is  simple,  i.e.  impulse,  step,  ramp,  or  sinusoid.  This  method  allowed 
manipulation  of  the  parameters  algebraically  in  the  frequency  domain,  instead  of  dealing  with 
differential  and  integral  Calculus  [24], 

The  other  method  of  solving  this  system  of  differential  equations  is  to  work  in  the 
frequency  or  ‘s’  domain  as  it  is  commonly  referred  to.  The  attractiveness  of  this  method  is  that 
the  solutions  can  be  obtained  using  algebra  to  solve  the  roots  of  the  characteristic  polynomial, 
which  represents  the  zeros,  poles,  and  gain  of  the  system  (ZPK).  Once  the  general  form  to  the 
solution  is  found,  eigenvectors  and  their  corresponding  eigenvalues  along  with  matrix  algebra 
and  a  table  of  Laplace  transforms  will  provide  a  solution  to  this  system.  Analysis  in  the  time 
domain  and  frequency  should  provide  the  same  result,  since  we  are  talking  about  the  same 
physical  system  under  analysis.  MATLAB™  is  equally  capable  of  working  in  this  domain  as 
well,  and  the  ‘tf2ss’  function  can  be  used  to  get  back  to  the  state-space  representation  if  desired. 

H(s)  is  the  transfer  function  and  maps  the  input  F(s)  to  the  output  Y(s)  in  the  frequency 
domain.  In  the  case  of  the  system  given  by  Eq.(27),  if  the  input  F  (t)  is  a  unit  impulse,  than  its 
Laplace  Transform  is  one  in  the  frequency  domain.  The  transfer  function  will  have  the  form  of: 

D 

1  sK  "thruster 

H(s)=  — - 1 -  (30) 

sz+2Cco  s+co  2 
n  n 

This  equation  will  generate  the  same  response  as  the  time  domain  state-space  representation  in 
equation  (27)  for  a  unit  impulse.  However,  since  all  the  data  and  parameters  are  measured 
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relative  to  time,  the  analysis  was  done  in  the  time  domain.  The  other  down  side  to  using  Laplace 
Transforms  to  solve  the  ODE  for  the  system,  is  that  the  inverse  Laplace  Transform  for  that 
forcing  function  must  be  known,  and  is  not  always  the  case.  The  first  situation  examined  was  the 
impulse  response  of  the  test  data  against  the  impulse  model  of  the  system. 

Test  Data  and  Models 


The  impulse  model  for  the  system  was  the  first  system  analyzed  to  see  the  system 
response  and  to  determine  how  closely  it  agree  with  data  for  a  measurement.  The  impulse  was 
supplied  by  the  firing  of  a  single  thruster  stick,  accomplished  by  rapidly  turning  the  thruster  off 
immediately  after  one  discharge.  The  state-space  of  the  system  is  found  by  writing  the  2nd  order 
ODE  as  a  set  of  first  order  ODEs  using  substitution  letting  Xi=0  and  X2=0,  the  system  has  the 
follow  representation  in  terms  of  matrices  A,  B,  C,  and  D.  A  is  called  the  state  transition  matrix, 
[24]  because  it  maps  the  states  from  one  value  to  another  as  the  input  signal  is  applied  to  the 
system 


0  A= 

x  =Ax+Bu,x=  ,u=F(t),  where  _-con2  -2qcon 

ej 

C=[l  0] 

y=Cx+Du 


B=  R. 


The  matrices  also  have  significance  in  the  frequency  domain  as  well,  since  the  transfer  H(s)  was 
computed  from  A,  B,  C,  and  D  through  the  formula: 


H(S)  =  (si  -  A) 1  •  B  +  D 


(32) 
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The  inverse  term  denotes  the  matrix  inverse.  The  above  system  can  also  be  built  in  Simulink™ 
which  was  used  to  verify  the  code  that  found  the  thrust  values,  and  to  check  out  the  system 
response  to  various  types  of  forcing  functions. 


This  is  tte  0.5Hi/l  47J  Model 


MOI.tolsJ- 01764 

Figure  3 1 .  Simulink™  model  for  Torsion  Balance  System  and  Thruster 


For  the  multi -pulse  model,  the  pulse  generator  was  set  to  operate  at  the  firing  frequency 
of  the  thruster  according  to  the  time  between  firings  provided  by  the  Cubloc™  software.  Both  in 
Simulink™  and  MATLAB™,  the  signals  usually  have  unit  amplitude.  This  effect  is  the  reason 
for  the  slide  gain  in  which  it  amplifies  or  attenuates  a  baseline  input  signal.  The  signal  is  then 
fed  into  the  state-space  model  to  obtain  the  output,  which  is  an  angular  deflection. 

The  impulse  response  for  1.47J  is  shown  in  Figure  32.  The  data  was  compared  against 
the  model  and  the  system  was  iterated  until  a  minimum  Root  Mean  Square  Error,  RMSE  of 
12.6197  grad  was  obtained.  This  corresponded  to  a  Ibit  of  approximately  107  pN-s  with  a 
standard  deviation  of  23  pN-s.  The  LDS  does  drift  a  little  over  time,  but  is  on  the  order  of  100 
pV,  and  does  not  significantly  affect  the  fitting  of  the  data  slice  to  the  model.  The  drift  is 
progressing  at  a  constant  rate,  so  it  can  be  removed  through  algebraic  techniques.  This  drift 
determination  essentially  involved  determining  the  linear  regression  coefficient  of  the  data  that 
had  a  drift  rate,  and  shift  the  data  in  the  opposite  direction  to  counteract  the  drift.  For  example,  if 
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the  data  had  a  slight  negative  slope  to  it,  the  data  would  be  shifted  by  a  positive  slope  of  equal 
magnitude  to  the  negative  slope  vectorially,  cancelling  out  the  drift  due  to  the  external  factors  or 
improper  C.G.  placement  or  external  vibration  . 


The  initial  set  of  data  was  taken  at  the  default  energy  discharge  setting  of  1.47J,  and  the 
nominal  operating  voltage  was  chosen  at  15  V,  the  midpoint  between  the  10-20V  ranges  stated  in 
operations  manual.  From  this  point,  the  voltage  varied  from  12.5  to  20  volts  in  increments  of  2.5 
volts.  The  original  plan  was  to  start  from  10  volts,  but  as  described  later,  10  volts  was  too  low  to 
operate  the  thruster  sticks  of  the  GG-pPPT  module. 

The  data  processing  required  the  use  of  a  Butterworth  filter  that  has  a  slope  of  -20*n 
decibels  (dB)  per  decade  (power  of  10),  where  n  is  the  order  of  the  filter.  A  fifth  or  sixth  degree 
filter  was  used  and  shown  in  Figure  33  and  is  used  as  low-pass  filter  for  the  data. 
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Gain=f{w/wc) 
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Figure  33.  Fifth-order  Butterworth  filter 

The  EMI  spikes  generated  by  the  discharge  capacitor  proved  to  be  very  useful  as  an 
accurate  clock  and  counter  of  pulses  from  the  thruster  sticks.  For  example,  using  the  spikes  from 
the  raw  data,  before  it  was  filtered  the  starting  point  for  the  multiple  pulse  data  set  was 
determined.  Because  of  the  very  short  nature  of  these  pulses,  the  pulses  would  accurately  depict 
where  or  when  the  thruster  fires.  This  ability  eliminated  the  need  to  rely  on  a  separate  detection 
sensor  to  count  the  pulses  as  the  thruster  fired  in  the  chamber.  The  voltage  spike  seen  from  the 
oscilloscope  were  much  greater  than  the  noise  floor  as  seen  in  (a)  of  Figure  34.  The  data  was 
further  processed  to  see  in  more  detail,  the  magnitude  and  response  of  the  balance  to  the  thruster 
as  shown  in  (b)  of  Figure  34. 
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(a)  EMI  spike  provides  a  reliable  way  to  determine  exact  firing  of  thruster,  within  the  resolution 

of  the  oscilloscope  and  sampled  data 


-3 


x  10 


(b)  Butterworth  filter  of  response  in  (a)  showing  the  detailed  oscillations  of  the  balance  from  the 

thruster  firing  at  the  set  frequency 

Figure  34.  Filtered  data  of  raw  system  response 

To  model  system  response,  a  model  of  the  input  forcing  function  requires  simulation, 

along  with  the  initial  conditions  of  the  system.  The  system  can  be  modeled  as  saw  tooth  or 
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square  wave  input.  The  type  of  forcing  function  model  used  will  depend  on  how  accurately  it 
represents  the  actual  shape  of  the  forcing  function  in  the  experiment. 

The  forcing  function,  F(t)  can  be  modeled  using  the  ‘gensig’  function  in  MATLAB™  or 
the  pulse  generator  in  Simulink™.  The  'gensig'  function  will  produce  a  unit  amplitude  pulse 
train.  It  must  be  multiplied  by  the  appropriate  gain  to  scale  the  forcing  function  to  the  correct 
value  needed  to  model  the  data.  The  multi-pulse  model  is  more  complicated  to  model  than  a 
single  impulse  model,  because  the  pulses  will  interact  with  each  other  preventing  a  complete 
decay  of  the  transient  condition  to  the  final  state  as  will  happen  with  an  impulse  input.  However, 
if  the  general  shape  of  the  forcing  function  is  known,  i.e.,  it  is  a  square,  sinusoidal,  exponential, 
linear  (ramp)  input  or  constant  (step)  input,  the  system  can  be  iteratively  solved  to  determine  the 
thrust  values  responsible  for  producing  those  waveforms. 

For  the  pulse  generator,  the  pulse  width  was  chosen  1.0  second  in  duration,  which 
allowed  the  force  determined  to  be  equal  to  the  average  Ibit  over  that  period  of  time.  The  actual 
pulse  trigger  discharge  occurs  very  quickly  and  based  on  discussions  with  Busek,  is  around  lOps. 
Thrust  generation  occurs  over  a  longer  time  frame  is  usually  in  the  range  of  60  to  80ps,  based  on 
plume  analysis  and  high-speed  imagery  captured  by  Selstrom  [13]  as  well  as  pulse  durations 
captured  during  this  research  effort.  Shown  in  Figure  35,  are  five  pulses  plotted  on  top  of  each 
other.  The  pulse  duration  in  Figure  35  was  determined  by  looking  at  the  point  where  the  voltage 
increased  above  the  baseline  noise  and  continued  upward  to  its  maximum  value,  and  then 
decreased  back  down  to  the  baseline  voltage  level,  which  was  4.03V  in  Figure  35.  This  pulse 
duration  represents  the  entire  ablation  an  ionization  timeframe.  The  actual  duration  of  the 
discharge  is  smaller  and  approximately  10  ps.  The  pulses  all  have  a  triangular  shape,  but  vary  in 
magnitude  and  duration  slightly  over  the  five  pulses  collected  and  plotted  in  Figure  35. 
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Pulse  1  - Pulse  2  — “Pulse  3  -Pulse  4  Pulse  5 


Figure  35.  PPT  pulse  shape  and  duration 

This  was  then  chosen  as  the  timeframe  initially  for  the  pulse  model  and  adjusted  up  to  1.0  second 
with  a  period  equal  to  the  inverse  of  the  firing  frequency  chosen.  Now,  we  will  look  at  some  of 
the  plots  generated  from  the  data  and  model  and  plot  them  on  top  of  each  other  for  comparison. 

The  1.47J  data  was  the  first  round  data  taken,  since  it  did  not  require  removal  from  the 
vacuum  chamber  and  reprogramming  of  the  thruster.  The  match  of  1 .0  Hz  data  with  the  model, 
shown  in  Figure  36  and  had  a  total  RMSE  of  8 1  grad  for  the  pulse  train.  The  data  under  ideal 
conditions  should  be  closer  to  the  shape  of  the  model,  but  even  the  use  of  the  Butterworth  filter 
was  not  able  to  remove  all  the  sources  of  noise  present  in  the  signal  or  perturbations  from  the 
nitrogen  pumps  that  would  turn  off  and  on  in  the  lab  during  testing.  Two  large  tanks  of  carbon 
dioxide  were  also  present  in  the  lab,  and  would  vent  from  time  to  time  to  regulate  the  pressure  in 
the  tanks.  Unfortunately,  when  this  occurred,  the  vibration  would  transit  from  the  tanks  to  the 
balance  through  the  floor,  affecting  the  measurements  slightly. 
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Figure  36.  1Hz  1.47J  data  and  model  plotted  on  same  graph 
The  Ibit  corresponding  to  this  plot  is  95.98  pN-s.  The  average  Ibit  over  all  runs  is  108 
pN-s.  Other  measurements  were  taken  and  the  thrust  and  Ibit  determined  from  solving  the  EOM 
for  the  system  and  the  forcing  function  respectively.  All  the  thrust  values  obtained  for  the  1.47J 
setting  for  various  frequencies  and  input  settings  are  shown  in  Table  5.  The  range  of  thrust 
values  tend  to  be  spread  out  for  this  sample,  but  are  higher  in  magnitude  than  the  rest  as  we  will 
see. 


Table  5.  Ibit  at  various  frequencies  and  voltages  for  1.47J 


Voltage 

Frequency  (Hz) 

Ibit  (pN-s) 

17.5 

2.00 

134.12 

15 

2.00 

95.98 

20 

1.50 

157.19 

15 

1.50 

105.99 

15 

1.00 

106.90 

12.5 

1.00 

95.98 

12.5 

0.67 

79.00 

17.5 

0.50 

107.48 

12.5 

0.50 

96.51 

Mean  Ibit 

(pN-s) 

108.79 

Standard 

Deviation 

(pN-s) 

23.35 
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The  next  set  of  plots  is  the  1.697J  data  collected  at  various  operating  voltages  and 
frequencies.  The  1.697J  setting  was  the  highest  programmable  setting  on  the  thruster,  for 
discharge  energy.  The  software  did  not  allow  the  selection  of  a  voltage  past  1300V,  without  risk 
of  damaging  the  unit,  and  a  two-pF  capacitor  provides  approximately  1.7J  of  energy  to  each 
thruster  stick  per  discharge  at  1300V.  The  match  between  the  data  and  the  model  was  iterated 
until  reaching  the  minimum  RMSE  was  found.  The  12.5V  1.697J  plot  is  the  first  one  in  this  set 
of  data.  The  RMSE  associated  with  the  run  is  12. 1  grad  and  is  shown  in  Figure  37.  The  data 
matches  the  model  well  in  this  case,  as  indicated  by  the  low  RMSE  value  compared  with  the 
previous  plot. 


From  (a)  of  Figure  37,  there  was  close  agreement  between  the  model  and  the  data,  the 
model  has  the  rougher  edges  because  of  the  sampling  period  of  the  simulation.  Smaller  sampling 
will  smooth  out  the  irregular  tips  of  the  model  plot.  The  first  50  seconds  are  shown  in  (b)  of 
Figure  37. 
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(b)  Close  up  of  the  first  50  seconds  of  the  pulse  train 


Figure  37.  1.697J  result  for  15V 
Table  6.  Ibit  at  all  voltages  and  frequencies  for  1.697J 


Voltage 

Frequency  (Hz) 

Ibit  (pN-s) 

17.5 

2.0 

105.68 

15.0 

2.0 

109.11 

12.5 

2.0 

70.11 

20.0 

1.5 

65.96 

17.5 

1.5 

71.17 

15.0 

1.5 

72.13 

12.5 

1.5 

109.23 

20.0 

1.0 

67.16 

17.5 

1.0 

77.07 

15.0 

1.0 

68.16 

12.5 

1.0 

60.73 

20.0 

0.5 

95.17 

20.0 

0.5 

66.15 

17.5 

0.5 

61.07 

17.5 

0.5 

54.83 

15.0 

0.5 

73.03 

15.0 

0.5 

62.24 

12.5 

0.5 

70.04 

12.5 

0.5 

52.30 

Mean  Ibit 

(pN-s) 

74.28 

Standard 

Deviation 

(pN-s) 

17.49 
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The  average  Ibit  for  the  1.697J  setting  is  74.28  pN-s  for  all  runs.  Looking  at  Table  6,  the 
Ibit  would  appear  to  increase  as  the  frequency  of  operation  of  the  thruster  unit  increases.  The 
reason  is  if  the  pulses  are  faster,  there  is  less  time  between  firing  and  the  propellant  face  does  not 
cool  as  much  as  it  would  at  the  lower  operating  frequencies.  Less  energy  is  required  to  ablate  the 
propellant,  still  hot  from  the  previous  pulse,  and  more  energy  will  go  into  ionizing  and 
accelerating  the  propellant,  to  a  high  Isp,  resulting  in  a  larger  Ibit  though  the  momentum  transfer 
of  the  ionized  Teflon™.  The  next  set  of  data  examined  is  the  1.6J  data. 


The  next  set  of  data  examined  was  the  1.6J  data  from  12.5  to  20V  and  shown  in  (a)  and 
(b)  of  Figure  38.  In  this  data,  the  ringing  and  oscillatory  nature  of  the  torsion  balance  for  an 
impulsive  thruster  is  evident 


(a)  1 ,6J  Entire  pulse  train,  showing  data  and  model 
The  oscillatory  ringing  observed  in  the  1 ,6J  pulse  train,  results  from  the  thruster 
providing  impulse  bits  in  resonance  with  the  motion  of  the  balance  and  increasing  until  the 
maximum  deflection  is  reached.  The  large  ‘hills’  in  the  data  show  the  effect.  The  average  Ibit 
over  all  frequencies  and  discharge  energies  is  36.1  pN-s  with  a  RMSE  of  12.51  prad.  The 
statistics  are  summarized  below  in  the  table  for  the  all  the  runs  at  1.6J. 
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(b)  1.6J  Pulse  train,  first  40s  showing  data  vs.  model  average  Ibit  is  36  pN 


Figure  38.  1.6J  Pulse  train,  comparing  model,  and  data 
Table  7.  Ibit  values  for  1.6J  at  all  voltages  and  frequencies 


Voltage 

Frequency(Hz) 

Ibit  (pN-s) 

12.5 

2.000 

55.35 

15.0 

2.000 

43.47 

17.5 

2.000 

42.96 

20.0 

2.000 

42.87 

12.5 

1.500 

29.73 

15.0 

1.500 

29.90 

17.5 

1.500 

16.19 

20.0 

1.500 

30.59 

12.5 

1.000 

36.29 

15.0 

1.000 

32.68 

17.5 

1.000 

35.71 

20.0 

1.000 

38.85 

20.0 

0.627 

38.19 

12.5 

0.500 

31.22 

15.0 

0.500 

34.71 

17.5 

0.500 

36.81 

20.0 

0.500 

37.76 

Mean  Ibit 

(pN-s) 

36.08 

Standard 

Deviation  (pN-s) 

8.20 
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From  this  table,  the  Ibit  does  have  a  slightly  higher  at  the  higher  voltages,  with  the 
exception  of  the  very  first  measurement.  The  one  unusual  thing  about  this  data  is  the  Ibit  was 
significantly  lower  than  the  thrust  observed  and  calculated  from  the  other  discharge  energies. 

The  exact  reason  is  unknown,  but  one  reason  could  be  that  the  sizing  of  the  pPPT  sticks  are  not 
favorable  for  this  discharge  setting  or  affected  by  the  time  constant  1/RC  of  the  circuit  powering 
the  thruster  sticks.  The  complete  spreadsheet  of  test  conditions  was  placed  in  appendix  B.  The 
next  set  of  data  will  consist  of  the  1.4J  data.  The  1.5Hz  data  is  in  Figure  39.  The  next  set  of  data 
that  was  collected  and  analyzed  is  the  1 ,4J  data,  and  a  plot  of  which  is  shown  in  Figure  39. 
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Figure  39.  1.5Hz  1.4J  at  20V 

The  RMSE  was  approximately  3.6  grad  for  the  1.4J  data  slice  in  Figure  39.  The  average 
Ibit  for  all  the  1.4J  0.5Hz  runs  was  approximately  92  pN-s,  and  the  average  over  all  runs  is  79.91 
pN-s  as  shown  in  Table  8.  This  result  is  consistent  based  on  the  period  of  the  calibration  pulse  of 
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1.357  seconds,  and  a  natural  frequency  of  4.63  rad/s.  The  deflection  is  lower  than  some  of  the 
other  plots  with  smaller  Ibit  values,  because  the  effective  spring  constant  is  so  much  larger  in  this 
case  requiring  a  greater  impulsive  force  for  the  same  amount  of  displacement  via  Hook’s  law  for 
a  torsional  spring  (F(t)*r=kspring0(t)).  The  interesting  aspect  of  this  data  is  the  fact  that  the  firing 
frequency  is  fast  enough  to  keep  the  arm  from  returning  fully  to  neutral  as  annotated  in  Figure  39 
Table  8.  Ibit  for  1.5Hz,  1.4J  data  at  all  voltages  and  frequencies 


Voltage 

Frequency 

(Hz) 

Ibit  (pN-s) 

12.5 

2.0 

49.70 

15.0 

2.0 

85.46 

17.5 

2.0 

88.60 

20.0 

2.0 

99.06 

12.5 

1.5 

91.80 

15.0 

1.5 

58.99 

17.5 

1.5 

55.20 

20.0 

1.5 

77.65 

12.5 

1.0 

96.53 

15.0 

1.0 

97.08 

17.5 

1.0 

87.87 

20.0 

1.0 

82.04 

12.5 

0.6 

33.83 

15.0 

0.6 

30.49 

17.5 

0.6 

63.30 

20.0 

0.6 

70.88 

12.5 

0.5 

107.89 

15.0 

0.5 

106.83 

17.5 

0.5 

107.52 

20.0 

0.5 

107.52 

Mean  Ibit 

(pN-s) 

79.91 

Standard 

Deviation 

(pN-s) 

24.25 
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The  other  phenomenon  observed  was  the  resonance  due  to  the  forcing  function  having 
the  same  frequency,  roughly  as  the  torsion  balance.  Previous  research  using  the  balance 
indicates  its  natural  frequency  should  be  in  the  0.62-0.67  Hz  range.  Resonance  was  very 
pronounced  at  two  frequencies,  one  at  0.668Hz  and  another  mode  present  around  0.627  Hz.  A 
graph  showing  the  phenomenon  of  the  0.668Hz  resonance  is  in  Figure  40.  In  this  graph,  the 
resonance  builds  to  the  point  where  it  is  almost  as  large  as  the  calibration  pulse  caused  by  the 
electrodes.  The  force  between  these  electrodes  is  919.22  pN  when  they  are  exactly  1.0  mm 
apart. 


Figure  40.  Resonance  at  0.668  Hz  when  firing  frequency  of  thruster  is  equal  to  natural  frequency 

of  torsion  balance 

For  the  resonance  to  build  to  this  magnitude,  the  forcing  function  of  the  thruster  is  exactly 
in  phase  with  the  oscillatory  motion  of  the  balance,  and  is  analogous  to  the  construct  interference 
and  resonance  of  wave  of  water  or  light  in  physics.  The  frequency  of  0.668  Hz  corresponds  to  a 
period  of  about  1.497s.  The  reason  for  wanting  to  know  the  resonance  frequency  is  that  the 
frequencies  and  modes  of  oscillations  are  typically  bad  for  structural  and  mechanical  stress 
reasons.  In  this  case,  however,  they  can  mask  the  true  effect  due  strictly  to  the  thruster 
operating.  Selstrom  [13]  estimated  that  the  natural  frequency  of  the  torsion  balance  should  be 
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around  0.66  Hz,  which  is  close  to  0.668.  The  balance  exhibited  some  modal  responses  similar  to 
Figure  40  but  did  not  reach  the  same  magnitude;  the  frequency  of  this  response  was  around  0.627 
Hz.  This  frequency  did  not  pose  any  problem  for  data  acquisition,  because  the  test  frequencies 
were  at  0.5,  1,  1.5,  and  2  Hz.  These  frequencies  are  far  enough  from  the  0.668  Hz  to  avoid  the 
excessive  oscillations  leading  to  resonance. 

The  next  set  of  data  to  analyze  was  the  1.3J  case.  Sometimes,  the  graph  needed  shifting 
left  or  right,  so  the  simulation  and  the  data  had  the  same  starting  point  at  t  equal  to  zero.  This 
shift  was  only  for  comparison  purposes  and  did  not  change  the  relationship  between  the  physical 
responses  of  the  torsion  balance  to  the  thruster  as  shown  in  Figure  41.  The  RMSE  was  3.37  grad 
for  this  sample  set.  With  some  of  the  graphs,  not  all  the  points  match  up,  because  of 
experimental  error  associated  with  the  equipment  and  taking  measurement. 

For  instance,  at  100  pV,  which  is  the  lowest  value  the  noise  level  reduced  to,  correspond 
to  an  angle  of  0.827  prad,  so  the  data  and  model  are  at  least  off  by  this  amount  at  a  minimum 
amount,  but  is  at  least  6  times  lower  than  the  value  typically  seen,  from  4.0  to  10  prad.  The 
average  Ibit  value  for  the  1.3J  setting  was  around  71  pN-s.  The  standard  deviation  for  this  data 
set  was  17.82  pN-s.  The  results  are  in  Table  9.  The  model  matched  the  data  closely,  once  the 
initial  conditions  are  considered,  and  added  to  the  simulation  as  an  initial  state  vector,  contain  the 
initial  angle,  0  and  the  angular  velocity,  0.  The  transient  response  of  the  model  compared  to  the 
data  was  very  sensitive  to  the  initial  conditions.  If  they  were  not  matched  properly,  then  the  two 
would  become  out  of  phase  much  quicker  than  anticipated  and  may  not  match  up  as  well,  since 
the  forcing  function  from  the  thruster  takes  time  to  imparts  its  frequency  characteristics  on  the 
torsion  balance  system. 
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Figure  41.  1.3J  data  at  0.5Hz  and  12.5V 
Table  9.  Ibit  for  1.3J  over  all  voltages  and  frequencies 


Voltage 

Frequency 

Ibit  (|iN-s) 

12.5 

2.0 

78.87 

15.0 

2.0 

84.16 

17.5 

2.0 

60.82 

20.0 

2.0 

42.87 

12.5 

1.5 

55.69 

15.0 

1.5 

96.21 

17.5 

1.5 

73.93 

20.0 

1.5 

82.49 

12.5 

1.0 

62.65 

15.0 

1.0 

39.97 

17.5 

1.0 

67.14 

20.0 

1.0 

82.43 

12.5 

0.5 

56.32 

15.0 

0.5 

95.27 

17.5 

0.5 

88.99 

Mean  Ibit 

OaN-s)  71.19 

Standard 

Deviation 

(jiN-s)  17.82 
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The  last  set  of  data  in  the  top  portion  of  the  data  matrix  is  the  1 .2  J  data  which  and  shown 
in  Figure  42.  The  Ibit  required  for  convergence  is  in  the  59  pN-s  range.  This  case  is  one  where 
the  drift  was  a  significant  hindrance  to  calculating  the  thrust.  The  drift  likely  caused  by  center  of 
gravity  error,  resulting  from  poor  placement  back  in  its  holder  on  the  base.  Fortunately,  it  is  a 
linear  function  of  distance  and  eliminated  mathematically  from  the  response  of  the  torsion 
balance  as  was  done  with  the  case  here.  This  setting  had  the  lowest  average  thrust  with  the 
exception  of  the  1.6J  data  shown  previously  in  Table  7.  The  results  of  the  1.2J  setting  are  shown 
Table  10. 


Table  10.  1 .2J  Ibit  at  all  voltages  and  frequencies 


Voltage 

Frequency 

Ibit  (pN-s) 

12.5 

2.00 

70.18 

15.0 

2.00 

91.72 

17.5 

2.00 

47.02 

20.0 

2.00 

92.25 

12.5 

1.50 

68.47 

15.0 

1.50 

38.72 

17.5 

1.50 

102.45 

20.0 

1.50 

62.64 

12.5 

1.00 

58.02 

15.0 

1.00 

42.86 

17.5 

1.00 

32.02 

20.0 

1.00 

57.39 

12.5 

0.50 

62.09 

15.0 

0.50 

41.10 

17.5 

0.50 

27.43 

20.0 

0.50 

58.70 

Mean  Thrust 

(pN-s) 

Standard  Deviation 

59.57 

(pN-s) 

21.85 
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Figure  42.  1.2J,  1.531Hz  15V 

Now  that  the  solutions  to  the  response  of  the  torsion  balance  can  be  accurately  modeled 
using  ODE  solution  techniques,  and  the  average  thrust  for  each  discharge  energy  setting  was 
found,  the  next  step  was  to  correlate  the  thrust  to  other  parameters  of  interest  like  the  discharge 
energy  frequency,  and  supplied  power.  We  will  now  make  some  observations  about  the  overall 
testing. 

All  the  Ibit  values  were  plotted  as  function  of  frequency  to  see  if  there  was  any 
correlation  between  the  Ibit  and  the  frequency  of  operation.  The  Ibit  values  seem  to  increase  in 
general  with  high  frequency,  but  more  samples  would  need  to  be  taken  at  the  various  frequencies 
to  verify  this  result.  A  scatter  plot  is  shown  in  Figure  43.  Statistically  speaking,  they  are  all  very 
close  to  the  same  value  and  have  similar  standard  deviations,  which  places  the  average  Ibit 
around  72  pN-s  for  all  thrust  values,  considering  discharge  energy  and  voltage  input.  The 
standard  deviation  for  the  each  data  is  larger  for  a  single  frequency,  than  between  frequencies. 
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Figure  43.  Chart  of  Ibit  values  vs.  frequencies 

The  other  objective  besides  examining  the  performance  of  the  thruster  due  to  variations 
in  discharge  energy  and  input  voltage,  was  to  attempt  to  determine  the  torques  produced  by  the 
sticks  separated  by  some  distance. 

Torque  Calculations 

The  other  goal  that  was  set  out  was  to  determine  the  torque  difference  if  any,  that  could 
be  measured  by  the  torsion  balance.  Due  to  time  restrictions,  only  data  at  1.3  and  1.2  J  were 
completed  out  of  the  original  planned  matrix. 

The  torque  was  calculated  using  the  difference  in  deflections  caused  when  stick  8  or  stick 
3  of  the  GG-pPPT.  Stick  8  was  the  stick  further  of  the  two  from  the  rotation  axis  of  the  balance. 
The  distance  between  the  stick  centers  is  0.585  in.  (.0149  m).  This  difference  in  length  is  what  is 
responsible  for  the  difference  in  torque  produced  between  the  two  sticks  assuming  that  the  thrust 
produced  is  the  same  for  both  sticks,  which  considering  the  close  agreement  of  earlier  data  is  a 
safe  assumption  to  make.  The  thrust  was  obtained  by  dividing  the  Ibit  by  1.0  s,  the  pulse  width 
used  in  the  model.  Neither  of  these  sticks  were  contaminated  with  carbon  build-up,  further 
increasing  confidence  in  this  assumption. 
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Table  11.  Torque  difference  produced  between  stick  8  and  3  at  1.3  J 


1.3J 

stick 

Voltage 

Thrust 

(pN) 

Rstick  (8)-m 

Rstick  (3)-m 

Torque(pN-m) 

Delta 

Torque 

(uN-m) 

8 

12.5 

92.20 

0.203 

0.191 

18.70 

3 

12.5 

81.50 

0.203 

0.191 

15.60 

3.127 

8 

15.0 

79.50 

0.203 

0.191 

16.10 

3 

15.0 

45.80 

0.203 

0.191 

8.74 

7.367 

8 

17.5 

59.60 

0.203 

0.191 

12.10 

3 

17.5 

35.60 

0.203 

0.191 

6.79 

5.284 

8 

20.0 

21.30 

0.203 

0.191 

4.32 

3 

Delta 

Torque 

(Avg)- 

pN-m 

20.0 

20.20 

0.203 

0.191 

4.09 

0.236 

4.003 

Table  12.  Torque  difference  produced  between  stick  8  and  3  at  1.2J 


1.2J 

stick 

Voltage 

Thrust 

Rstick  (8)-m 

Rstick  (3)-m 

Torque(pN-m) 

Delta 

Torque 

(uN-m) 

8 

12.5 

86.70 

0.203 

0.191 

17.60 

3 

12.5 

91.30 

0.203 

0.191 

17.40 

0.122 

8 

15.0 

80.40 

0.203 

0.191 

18.90 

3 

15.0 

99.20 

0.203 

0.191 

16.30 

2.638 

8 

17.5 

95.30 

0.203 

0.191 

20.20 

3 

17.5 

106.00 

0.203 

0.191 

19.30 

0.924 

8 

20.0 

71.00 

0.203 

0.191 

21.50 

3 

Delta 

Torque 

(Avg)- 

pN-m 

20.0 

75.40 

0.203 

0.191 

14.40 

7.08 

2.691 
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From  Table  1 1  and  Table  12,  there  is  a  small  difference  between  the  two  but  hardly 
noticeable  and  requires  the  noise  floor  and  external  vibrations  to  be  as  minimal  as  possible. 
Though  a  torque  produced,  it  is  small  compared  to  the  actual  magnitude  of  the  thrust  itself.  With 
the  average  for  1 .3  J  at  4  jiN-m  of  torque,  the  average  thrust  (Ibit/pulse  duration)  is  somewhere 
between  70  and  100  pN  and  amounts  to  only  six  percent  of  the  thrust  produced  by  each  pPPT 
thruster  stick.  These  thrust  levels  are  not  large  enough  to  rotate  a  CubeSat  with  a  mass  of  few 
kilograms.  It  is  safe  to  say  that  the  separation  between  the  sticks  on  the  GG-pPPT  module  will 
not  result  in  a  nutation  of  the  CubeSat. 

Observations 

There  were  several  important  things  discovered  about  the  operation  of  the  thruster  as  it 
was  firing.  The  first  thing  found  was  the  thruster  has  a  minimum  voltage  required  for  operation 
higher  than  the  lowest  value  stated  in  the  operations  manual. 

Discovery  number  one,  the  GG-pPPT  minimum  operation  voltage  is  between  10  and  12.5 
volts.  During,  preliminary  investigation  of  the  thruster,  the  input  voltage  was  turned  to  the 
lowest  recommended  setting  of  10V.  Everything  else  on  the  thruster  was  powered  on  and  ready 
for  operation.  When  the  frequency  knob  on  the  stick  selector  module  was  turned  on,  none  of 
sticks  fired.  Had  these  sticks  fired,  the  Cublock™  software  and  the  firmware  in  the  thruster 
would  have  relayed  a  signal  indicating  a  fire  or  misfire  of  the  thruster  as  shown  in  the  screen  in 
Figure  44.  This  screen  contained  other  useful  information  about  the  behavior  of  the  GG-pPPT 
thruster  module.  For  example,  it  would  display  the  time  between  firings,  essential  in 
determining  the  frequency  of  the  forcing  function  that  went  into  the  MATLAB™  and 
Simulink™  models.  It  also  contained  very  useful  information  like  the  discharge  voltage  of  the 
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main  capacitor  and  the  trigger  capacitor,  it  tracked  system  time  and  the  number  of  misfires  for  a 
given  pPPT  stick. 


Figure  44.  Screenshot  of  cublock™  debug  terminal 
In  this  screenshot  for  instance  the  voltage  of  the  trigger  capacitor  is  7845  volts,  and  the 
main  capacitor  is  discharging  at  1213  V,  about  1.47J  discharge  energy.  There  would  have  also 
been  a  flash  of  light  that  would  have  been  observed  through  the  side  viewport  on  the  left  side  of 
the  vacuum  chamber  as  well,  and  would  have  looked  like  bright  blue  flashes  shown  in  Figure  45 
This  result  is  important  for  the  satellite  this  device  is  intended  for,  because  of  voltage 
requirements.  The  exact  cutoff  voltage  for  operation  was  not  determined,  but  operation  at  12.5 
volts  allowed  continuous  thruster  module  operation.  So  the  limit  is  somewhere  between  10  and 
12.5  volts.  The  flashes  of  light  were  also  dimmer  at  the  lower  voltage  settings  of  12.5  volts  vs. 
20  volts.  The  viewing  stand  shown  in  the  equipment  section  in  (a)  of  Figure  27,  made  it 
relatively  easy  to  see  the  operation  of  the  thruster  in  the  confined  space  of  the  vacuum  chamber. 
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At  times,  there  were  instances  when  the  sticks  would  misfire  or  two  to  three  of  them  would  fire 
rapidly  in  succession.  The  thruster  side  window  was  on  the  left  side,  looking  towards  the  front 
of  the  chamber. 


(a)  Orientation  of  thmster/mirrors  (b)  High-resolution  of  picture 

Figure  45.  Side  profile  operation  of  GG-PPT 


Discovery  number  two  was  while  the  thruster  was  operating,  not  all  the  sticks 
accumulated  carbon  during  operation,  or  at  the  same  time  for  those  sticks  where  carbon  deposit 
buildup  did  occur.  The  first  stick  to  have  a  visible  carbon  buildup  on  the  Teflon™  surface  was 
stick  number  seven.  Each  time  the  thruster  was  removed  out  of  the  vacuum  chamber  for 
reprogramming,  pictures  of  the  thruster  and  thruster  heads  were  taken.  Each  thruster  head  was 
inspected  to  check  for  cracks  in  the  surface  or  damage  to  any  of  the  electrode  surfaces.  In 
general,  the  center  electrode  and  cathode  receded  away  with  the  Teflon™  as  expected.  The 
surface  however  did  not  remain  entirely  flat,  some  of  the  thruster  sticks  had  a  ‘bowl’  or  curved 
Teflon™  surface.  The  sticks  were  still  firing  reliably  throughout  testing. 

Looking  at  (b)  of  Figure  46,  the  bowl  or  conical  depression  in  the  Teflon  is  visible.  The 
inner  electrodes  visible  in  the  center  of  each  stick,  also  recedes  with  the  Teflon™  as  designed 
during  thruster  operation.  The  initial  thought  was  the  thruster  plume  could  be  returning  to  the 
surface  of  the  thruster,  because  it  is  in  close  proximity  to  the  wall  of  the  chamber.  If  that  was  the 
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only  source  then  all  the  thrusters  should  have  equal  amounts  of  carbon  on  them  statistically 
speaking. 


I  Carbon  buildup 
Stick  5  and7 


(a)  Close-up  of  stick  7  (b)  Contamination  now  present  on  sticks  5  and  7 

Figure  46.  Carbon  build-up  on  thruster  heads 
From  (a)  in  Figure  46,  the  carbon  build  up  on  this  sticks  is  shown  by  the  arrow, 
sometimes  the  thrusters  would  discharge  prematurely  or  more  rapidly  than  would  be  expected, 
based  on  the  set  firing  frequency  for  the  thruster  module.  These  carbon  deposits  are 
semiconductors  when  they  are  hot  and  vaporized  and  could  lead  to  some  shorting  of  the  potential 
between  the  cathode  and  the  inner  electrode.  The  other  noticeable  effect  on  the  data  was  if  the 
thruster  center  of  gravity,  C.G  was  either  forward  or  back  of  the  counterweight  C.G.  This  effect 
showed  up  as  a  positive  or  negative  drift  of  the  waveform  observed. 

Discovery  number  three,  the  high  drift  rates  are  due  to  CG  location  error  and  tension 
changes  in  wiring.  If  the  drift  rates  were  thermal  in  nature,  they  would  likely  have  averaged  out 
to  about  the  same  amount  each  time.  However,  some  of  the  rates  were  has  high  as  10  prad  over 
a  100  s  period  indicating  the  source  was  due  to  something  else.  The  drift  rates  were  also 
consistently  negative.  The  wires  had  some  associated  weight  with  them  and  not  easily  accounted 
for  in  the  center  of  mass  calculations  discussed  earlier. 
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Figure  47,  provides  an  example  of  the  drift  present  in  some  of  the  data  that  had  to  be 
corrected.  The  figure  depicts  the  raw  change  in  LDS  voltage  vs.  time,  and  the  negative  drift  due 
to  mirror  moving  further  from  the  LDS  tip  is  evident  from  the  negative  constant  slope  of  the  data 
in  Figure  47.  However,  because  it  was  linear,  it  was  a  simple  correction  factor  to  make  to  the 
original  data  file  and  turn  it  into  a  graph  similar  to  Figure  48.  This  technique  was  applied  to  the 
data  for  accuracy  and  consistency,  had  to  done  to  the  files  earlier  allow  the  plot  of  the  data 
against  the  model  to  make  sense.  Once  adjusted,  the  data  could  then  be  compared  directly 
against  the  ODE  model  for  the  system,  and  the  thrust  corresponding  to  the  minimum  RMSE  with 
the  same  timing  frequency  as  the  thruster  could  be  found.  Initial  thoughts  were  the  thruster  may 
be  forcing  the  drift.  If  the  thruster  was  the  cause,  the  drift  would  have  been  in  the  opposite 
(positive)  direction,  the  direction  the  the  torsion  arm  moves,  moving  the  mirror  closer  to  the  LDS 
tip.  Correcting  the  data  for  the  center  of  gravity  error  allows  the  MATLAB™  to  assume  that 
there  is  no  placement  error  when  running  the  simulation.  The  result  is  the  ability  for  the  code  to 
converge  and  find  a  minimum  RMSE  value  for  the  specific  thrust  accurately. 
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Figure  47.  Uncorrected  data  with  drift  rate 
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Figure  48.  Corrected  data  without  drift 

Discovery  number  four,  the  thrust  is  primarily  a  function  of  discharge  energy  and  firing 
rate.  The  voltage  did  seem  to  increase  slight  at  20  volts,  but  this  was  not  always  consistent 
across  discharge  energies.  In  some  cases,  for  example  the  12.5  volt  setting  produced  more  thrust 
than  the  20  volt  setting  for  the  same  discharge  energy.  However,  the  chart  in  Figure  43  indicates 
that  there  is  a  small  but  noticeable  increase  in  the  average  thrust  as  the  firing  rate  is  increased. 

This  is  increase  was  unexpected  initially.  The  power  supply  will  maintain  constant 
voltage,  and  the  capacitor  will  only  charge  to  the  voltage  set  by  the  discharge  energy  setting.  So, 
increasing  the  voltage  without  increasing  discharge  energy  has  little  effect  on  the  calculated 
thrust  from  the  thruster.  However,  if  the  voltage  is  insufficient,  then  the  electronics  cannot 
amplify  the  given  voltage  to  the  level  needed  for  discharge,  and  was  seen,  when  powered  at  only 
10V,  the  thrusters  would  not  fire  at  all. 


82 


Error  Analysis 

This  section  on  the  error  analysis,  identifies  the  various  sources  of  error  in  the  research 
experiment.  We  begin  by  discussing  the  noise  present  in  the  signal  coming  from  the  LDS,  and 
the  signal  to  noise  ratio.  From  there  the  error  associated  with  the  sample  time  and  time  steps 
chosen  are  discussed,  and  then  measurement  error  in  the  output  of  the  Agilent  54622D  Mixed 
Signal  Oscilloscope.  Next,  we  generate  some  normal  probability  and  histogram  plots  for  the 
thrust  data.  We  end  this  section  with  discussion  on  the  norm  of  the  residuals 
Signal  to  Noise  (SNR) 

The  biggest  source  of  noise  was  the  high  frequency  turbo-pump  and  the  surrounding  N2 
pumps  that  were  co-located  in  the  lab  where  the  research  was  performed.  The  signal-to-noise, 
SNR  ratio  is  a  useful  quantity  to  determine  how  much  stronger  the  signal  is  than  the  noise  floor 
of  the  environment.  Initially  the  oscilloscope  was  used  as  is.  Early  attempts  were  made  to  collect 
quality  data,  but  the  noise  floor  was  around  lOmV  and  much  too  high  to  accurately  see  the 
motion  of  the  balance  due  to  thruster  operation  without  extensive  signal  processing  techniques. 

Most  of  the  pulses  on  average  were  only  one  or  two  mV  in  magnitude.  Operating  key 
features  of  the  oscilloscope  like  the  high  frequency  noise  rejection  mode,  eliminated  sources 
greater  than  50Hz  in  frequency  greatly  reducing  the  noise  present  in  the  signal.  The  data  sample 
was  averaged  using  the  previous  point,  and  resulted  in  a  cleaner  signal,  and  eliminated  much  of 
the  random  noise  present. 

With  these  settings,  the  noise  floor  was  reduced  from  10  mV  to  approximately  100  to  260 
pV.  This  corresponds  to  an  increase  in  detection  sensitivity  of  38X.  This  was  especially 
important  in  this  research,  because  of  the  relatively  massive  thruster  and  counterweight.  Other 
individuals  had  thrusters,  but  one  of  them  was  only  one-tenth  the  mass  of  the  GG-pPPT  module. 
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An  increase  by  a  factor  of  38  is  easily  converted  to  decibels,  dB  and  is  the  common  unit  way  of 
measuring  the  power  and  intensity  of  signals  in  the  electrical/acoustical  engineering  world.  The 
value  in  dB  was  found  using  the  formula: 

dB=10  log^21!L  =  —  (33) 

P  I 

in  in 

P0, 10  is  the  power  or  intensity  output  and  P;n,  Im  is  the  power  or  intensity  input  respectively. 
Substituting  numbers  into  the  equation,  the  value  corresponding  to  this  change  was  from  10  mV 
to  260  pV  is  +15.8dB  for  example. 

Time  Step  Resolution 

The  impulsive  nature  of  the  thruster  and  the  relatively  short  pulse  of  the  thruster  and 
oscillations  required  instrumentation  capable  of  measure  brief  instants  in  time  but  enough 
memory  capacity  to  store  a  couple  of  minutes  of  data  during  a  experimental  run.  The  Agilent 
54622D  Oscilloscope  has  a  dynamic  range  from  5.0  ns  to  50.0  s.  This  range  covers  the  fastest  of 
familiar  events.  At  5.0  ns,  even  a  beam  of  light  travelling  at  670  million  mph  will  only  travel 
about  5  ft!  The  original  data  files  had  two  to  four  million  data  points,  but  was  discovered  to 
cause  problems  during  data  analysis,  because  of  the  time  required/and  computing  resources.  A 
decision  was  made  to  reduce  this  number  to  4000  data  points  per  run,  which  Excel  could  handle 
for  quick  assessment  on  the  quality  of  the  data. 

Measuring  Instrument  Error. 

The  error  associated  with  taking  the  masses  of  the  counterweight  are  +/-  0.1  g  for  the 
scales  that  were  used  in  ASYS  632-Satellite  Design  and  Test.  The  calipers  used  for  measuring 
the  distances  on  the  torsion  balance  had  accuracy  down  to  0.01  mm  or  0.397  mil.  The  Agilent 
54622D  oscilloscope  had  a  dynamic  range  in  voltage  of  5  V  to  1  mV/division. 
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Data  Error  and  Statistics 


The  statistics  have  been  tabulated  for  all  the  Ibit  data  and  shown  in  and 
Table  13.  The  maximum  Ibit  appeared  to  be  around  1.47J. 

Table  13.  Mean  and  standard  deviation  statistics 


Discharge 
Energy  (J) 

Mean  Ibit 
(uN-s) 

1.47 

108.79 

1.69 

77.43 

1.6 

35.54 

1.4 

81.19 

1.3 

71.56 

1.2 

59.57 

Total 

72.35 

Standard  Dev.,o 

(pN-s) 

1.47 

23.35 

1.69 

17.49 

1.6 

7.82 

1.4 

8.20 

1.3 

17.82 

1.2 

21.85 

Total 

16.09 

The  largest  standard  deviation  is  from  the  1.47J  data  with  a  standard  deviation  of  23.35 
pN-s.  All  other  values  have  standard  deviations  that  are  smaller.  The  total  average  Ibit  for  this 
thruster  over  all  energies  and  frequencies  is  approximately  73  pN-s  and  comparable  to  other 
micro-thrusters  with  similar  power  requirements.  Busek  tested  these  sticks  at  1.96J,  and 
indicated  a  average  Ibit  153  pN-s  as  the  starting  max  for  the  pPPT  sticks,  and  80  pN-s  was  the 
average  throughout  the  life  of  the  thruster  until  all  propellant  has  been  consumed  according  to 
the  operations  manual  and  shown  in  Figure  49.  The  average  mass  loss  per  pulse  was  around 
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19.71  pg/pulse,  according  to  the  manufacturer,  with  a  Isp  of  827s  [25].  The  Isp  corresponds  to  a 
velocity  of  8110.1  m/s  or  26,609  fps. 


Figure  49.  Manufacture's  Ibit  performance 
Table  14.  Summary  statistics  for  all  data 


Parameter 

Min 

Max 

Max  Ibit-pN-s 

157.00 

Ibit  avg. 

(all  data)-pN-s 

72.75 

Standard 

Deviation 

14.00 

24.00 

(avg./max)-pN-s 
Number  of 
firings 

Test  Conditions 

24216.00 

112.00 

Torque(pN-m) 

1.3J 

4.00 

Torque  (pN-m) 
1.2J 

2.69 

Power 

Consumed 

31.80 

35.40 

(pulsed 

min/max)-W 
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Thrust  efficiency  was  determined  once  the  discharge  energy,  pulse  time,  the  velocity  and 
mass  loss  rate  are  known.  For  example,  at  1.47J/shot  and  using  the  mass  loss  of  19.7 1*1  O’9 
kg/shot  and  ue  equal  to  81 10.1  m/s,  results  in  a  theoretical  efficiency  of  about  31  percent  using 
equation  (12).  This  is  just  the  thruster  efficiency.  The  total  efficiency  will  be  less  than  this 
because  of  the  all  the  heat  losses  of  the  electronics  in  inside  the  GG-pPPT  and  the  heat  lost  to 
producing  hot  neutral  particles,  not  accelerated  by  the  Lorentz  Force.  The  thruster  efficiency 
was  determined  to  be  about  21%  in  this  research,  but  that  does  not  include  the  losses  from  the 
internal  electronics  of  the  GG-pPPT  device.  Factoring  these  in,  the  actual  efficiency  is  probably 
close  to  historical  values  of  6-10  percent. 

In  summary,  the  errors  of  the  instrumentation  and  the  test  collection  tools  were  at  least 
four  to  10  times  smaller  than  the  signal  analyzed  and  collected  from  the  oscilloscope.  Moving 
the  roughing  pump  off  the  stand  and  using  the  noise  rejection  features  of  the  oscilloscope 
resulted  in  a  much  better  SNR,  making  data  analysis  and  model  fitting  much  easier.  The  high 
frequency  noise  is  random  and  in  nature,  and  removed  using  filtering  techniques.  The  error 
between  data  collected  in  this  research  and  Busek  is  2.6  %  for  the  maximum  value  and  9.1%  for 
the  average.  . 

V.  Conclusions  and  Recommendations 

The  final  chapter  of  this  research  take  the  information  gathered  from  the  thrust,  plots, 
photographs,  error  plots  and  statistical  fitting  to  determine  the  best  conditions  to  operate  at,  and 
the  test  conditions  that  were  accomplished.  Next,  recommendations  are  made  to  individuals 
following  up  with  this  research  or  upgrading  the  torsion  balance  to  be  adjustable  once  down  to 
vacuum  conditions.  Finally,  there  is  some  discussion  regarding  the  future  work  on  this  device. 
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Optimum  Voltage  and  Energy  Operating  Conditions 

The  first  major  finding  is  the  Ibit  appears  to  decrease  with  a  decrease  in  discharge  energy. 
However,  from  the  graphs  presented  in  chapter  four,  some  variation  with  the  data  as  is  present  as 
the  conditions  change  and  the  thruster  moves  in  and  out  of  the  chamber  during  reprogramming, 
and  effectively  change  the  spring  constant.  This  is  a  concern  because  spring  constant  that  is  too 
stiff  will  prevent  the  signal  from  being  above  the  noise  floor.  The  total  average  for  Ibit  across 
all  runs  is  73  pN,  close  to  spec.  The  best  setting  to  operate  the  GG-pPPT  is  20V  1.5Hz  at  1.47J 
for  maximum  Ibit  of  157  pN-s.  The  next  best  setting  is  at  1.4J  and  provides  approximately  83 
pN-s  of  impulse. 

Earlier  Ibit  levels  that  were  calculated  were  about  181  pN-s,  but  that  number  was  slightly 
too  high  for  two  reasons.  First,  the  sensitivity  constant  of  -0.4V/mm  or  -400  mV/mm  that  was 
used  is  slightly  too  low.  Using  the  chart  provided  by  Philtec™,  the  sensitivity  constant  of -15.30 
mV/mil  for  the  far  side  of  the  LDS  sensor  translates  into  -0.602V/mm  or  -602  mV/  mm.  This 
sensitivity  will  affect  the  observed  voltage  output  and  change  on  the  oscilloscope  and  will  make 
thrust  levels  appear  higher,  than  they  really  are. 

Second,  even  though  the  thrust  force  is  relatively  constant,  the  impulse  bit  over  time  will 
decrease  due  to  a  larger  proportion  of  particles  not  ionized  as  efficiently  as  when  the  thruster  was 
new.  The  thruster  efficiency  calculated  solely  on  the  discharge  energy  results  in  a  theoretical 
efficiency  of  3 1  percent.  The  main  culprit  to  reduced  efficiency  would  be  the  particle 
accumulation  that  accrues  on  the  surface  of  the  Teflon™.  This  effectively  reduces  the  Isp  and 
uses  more  propellant  resulting  also  in  decreased  efficiency. 
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Not  all  the  conditions  for  the  torque  were  accomplished,  but  torque  was  measured 
between  two  different  energy  levels  at  all  four  voltage  settings  for  the  1.3  and  1.2J  energy  setting 


in  Table  15. 


Table  15  Accomplished  Test  Conditions 


Control  variables 
Discharge 
energy  (J) 

E=l/2CV2,C=2gF 

Output  voltage 

Frequency(Hz) 

Stick  pattern 

Voltage 

S 

1.7 

1303.84 

0.5, 1,1. 5, 2 

1,2, 3, 4, 5, 6, 7, 8, 9 

12.5,15,17.5,20 

S 

1.6 

1264.91 

0. 5,1,1. 5, 2 

1,2,3, 4, 5,6, 7, 8,9 

12.5,15,17.5,20 

V 

1.5 

1224.74 

0.5, 1,1.5, 2 

1,2, 3, 4, 5, 6, 7, 8, 9 

12.5,15,17.5,20 

V 

1.4 

1183.22 

0. 5,1,1. 5, 2 

1,2,3, 4, 5,6, 7, 8,9 

12.5,15,17.5,20 

s 

1.3 

1140.18 

0.5, 1,1.5, 2 

1,2, 3, 4, 5, 6, 7, 8, 9 

12.5,15,17.5,20 

1.2 

1095.45 

0. 5,1,1. 5, 2 

1,2,3, 4, 5,6, 7, 8,9 

12.5,15,17.5,20 

Discharge 
energy  (J) 

Output  voltage 

Frequency(Hz) 

Stick  pattern 

Voltage 

1.7 

1303.84 

0.5, 1,1.5, 2 

3,8 

12.5,15,17.5,20 

1.6 

1264.91 

0. 5,1,1. 5, 2 

3,8 

12.5,15,17.5,20 

1.5 

1224.74 

0.5, 1,1.5, 2 

3,8 

12.5,15,17.5,20 

1.4 

1183.22 

0. 5,1,1. 5, 2 

3,8 

12.5,15,17.5,20 

S 

1.3 

1140.18 

0.5, 1,1.5, 2 

3,8 

12.5,15,17.5,20 

1.2 

1095.45 

0. 5,1,1. 5, 2 

3,8 

12.5,15,17.5,20 

Time  was  an  issue,  because  of  the  need  to  share  equipment,  but  there  is  no  reason  to 
suspect  that  the  torque  difference  between  the  sticks  would  be  any  different  for  the  other 
discharge  energy  levels,  because  the  relative  thrust  is  the  same  for  all  sticks,  unless  contaminated 
or  damaged.  The  firing  of  the  pPPT  sticks  was  dependable,  but  there  was  misfiring  that 
occasionally  did  occur.  One  possible  explanation  for  the  misfiring  is  that  the  plume  ion  plume 
from  the  adjacent  stick  is  sufficient  to  effectively  short  out  the  electric  potential  established  by 
electric  field  that  exists  between  the  various  electrodes  in  the  thruster  unit  and  the  cathode.  The 
cathode  is  rated  to  reach  10  kilovolts,  9,700  volts  was  the  highest  observed  in  research.  With 
such  a  high  electrical  potential,  even  a  good  vacuum  might  allow  enough  particles  to  start  an 
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ionization  trail  in  addition  to  the  charged  particles  coming  out  of  the  other  thruster  sticks,  at  the 
tight  spacing  of  this  assembly. 

The  very  fast  response  of  the  oscilloscope  allowed  the  capture  of  details  like  the  pulse 
width  of  the  thruster  firing  that  would  have  been  difficult  to  measure  otherwise.  There  are  a  few 
recommendations  put  forth  for  the  next  generation  of  torsion  balance  or  for  an  individual  doing 
follow-on  research  with  the  GG-pPPT. 

Recommendations 

Torsion  Balance  Improvements  and  Considerations 

There  is  much  to  say  about  the  torsion  balance  provided  by  Busek.  The  sensitivity  of  the 
instrument  is  amazing  to  say  the  least.  Very  few  devices  can  measure  such  small  impulse  bits 
and  forces  produced  by  pPPTs  and  other  'micro-thrusters'  accurately,  even  with  some 
background  noise  present.  There  are  a  couple  of  improvements  to  consider  for  future  models  or 
upgrades.  First,  I  would  start  with  the  calibration  process  for  the  main  and  damping  electrode. 
Under  ideal  conditions,  the  distance  between  the  electrodes  should  not  change  significantly  from 
the  1 .0  mm  position  if  not  disturbed.  However,  the  reality  of  the  situation  is  different.  The 
calibration  process  takes  place  in  the  external  environment  and  atmosphere,  and  thermal 
gradients  as  well  as  physical  changes  of  the  sized  due  to  pressure  forces,  change  this 
distance..  When  the  chamber  is  under  vacuum,  the  dimension  of  the  chamber  due  to  thermal  and 
atmospheric  effects  can  cause  the  alignment  of  the  balance  to  change,  where  the  electrodes  are 
no  longer  1  mm  apart.  This  will  change  the  voltage  to  force  calibration  constant. 

A  recommended  action  would  be  to  attach  a  separate  low  noise  motor  be  connected  to 
the  micrometer  with  electrical  connections  so  that  the  distance  of  1.0  mm  could  be  re-established 
once  in  vacuum  conditions  and  the  system  is  stable  in  terms  of  temperature  and  physical  strain. 
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The  other  change  recommended,  would  be  the  installation  of  some  electrical  connections 
on  the  base  of  the  balance  arm.  One  of  the  toughest  issues  was  dealing  with  the  26  or  so  wires 
that  were  very  stiff  and  overwhelmed  the  spring  constant  of  the  balance  itself.  This  would  reduce 
the  amount  of  error  in  both  the  calculation  and  the  data  gathering.  A  standard  serial  connection 
and  power  connections  would  have  been  immensely  helpful.  All  in  all,  this  balance  is  the  most 
sensitive  thrust  measuring  instruments  at  AFIT.  The  instrument  was  so  sensitive  that  most  of  the 
data  gathering  had  to  be  done  late  in  the  evening,  when  the  amount  of  external  noise  was  as  low 
as  it  could  be  (minus  the  N2  pumps.) 

Future  Research 

One  variable  for  further  investigation  is  the  duty  cycle  of  both  the  main  capacitor  and  the 
trigger  circuit  on  the  Ibit.  The  default  is  30  percent  for  the  trigger,  though  at  its  maximum 
energy  per  pulse  of  1.697  J,  the  main  capacitor’s  duty  cycle  jumped  to  100  percent.  There  was  no 
information  on  the  effect  of  this  parameter  on  the  overall  performance  of  the  thruster.  Changing 
this  parameter  without  damaging  the  internal  electronics  or  the  thrusters  themselves  was  a 
principle  concern,  since  this  is  a  one  a  kind  piece  of  equipment.  I  learned  a  great  deal  about 
electronics,  plume  interaction  and  EMI,  through  my  research  with  this  thruster  was  amazed  by 
how  brilliantly  it  shines.  This  is  definitely  a  propulsion  system  that  any  nano- satellite  program 
would  benefit  from. 

In  summary  the  highest  Ibit  observed  was  157  pN-s,  with  a  average  of  72  pN-s  across  all 
discharge  energies.  In  some  of  the  plots,  the  Ibit  seemed  to  be  slightly  higher  for  the  higher 
voltages,  and  in  others,  it  was  not  the  case.  Firing  frequency  had  a  more  dramatic  effect  on  the 
Ibit  than  the  voltage,  but  overall  the  Ibit  varies  the  most  by  discharge  energy  as  evident  in 
Table  13. 
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The  photo  in  Figure  50,  is  the  GG-pPPT  facing  the  window,  taken  with  a  SLR  camera 
with  a  exposure  duration  of  five  seconds,  because  at  two  hertz  firing  frequency,  all  thrusters  can 
be  captured  in  the  photo. 


Figure  50.  Gatling-gun  firing  at  2  Hz 
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Appendix  A.  Physical  Dimensions  of  the  Torsion  Balance 


Thesis:  Gatling  Gun-micro  Pulsed  Plasma  Thruster  (PPT) 

24  Sept  2010-Equipment  setup-calibration 
Thruster  Physical  Characteristics  and  Measurements 
Mass:  1 .898  Kg  (Satellite  Design  Scale) 

1.895  Kg  (Chemical  Back  Lab  Scale) 

Dimensions: 

(L/H)  7-7/32"~7.25  in  (Length  or  height) 

(W)  4-7/32"~4.25  in  (Width-side  to  side) 

(D)  5-7/32"~5.25  in  (front  to  back) 


Center  of  Mass  Setup: 


Center  of  Mass  properties 

Zbar  (L/H):  3-29/32"  from  top 

Ybar  (from  power  supp.  side)  2-6/32-2.1875" 

Xbar  (from  front  to  back 


Counter  weight 

•  Material:  T6061  A1 

•  V=45  in.3 

•  M=1.921  kg 

•  10/32  hole  drilled  base  of 
weight  for  mounting 


Torsional  Balance  Characteristics 

Arm  length :  1 8-24/32=1 8.75"  edge-edge 

Geometric  Center:  1 ,5"x  1 .5"  on  each  plate 

Distance  between 

geometric  centers:  1 5-24/32=1 5.75" 

Hole  spacing  (arm-arm) 

15.1225mm~0.5954" 

Hole  spacing  (side-side) 

8.6225mm~0.3395" 
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Appendix  B:  Test  Conditions 


Test 

Condition  # 

Frequency 

(Hz) 

Discharge  Energy 

(J) 

1. 

0.500 

1.47 

2. 

0.500 

1.47 

3. 

0.668 

1.47 

4. 

1.000 

1.47 

5. 

1.000 

1.47 

6. 

1.500 

1.47 

7. 

1.500 

1.47 

8. 

2.000 

1.47 

9. 

2.000 

1.47 

10. 

0.500 

1.697(No  Cal) 

11. 

0.500 

1.697 

12. 

0.500 

1.697 

13. 

0.500 

1.697 

14. 

0.500 

1.697  (No  Cal) 

15. 

1.697 

1.697  (No  Cal) 

io. 

Ui  j  uu 

17. 

0.500 

1.697(No  Cal) 

18. 

1.000 

1.697 

19. 

1.000 

1.697 

20. 

1.000 

1.697 

21. 

1.000 

1.697 

22. 

1.500 

1.697 

II 

ffil 


2.000 


0.500 


0.500 


25 

..5.. 

1.000 

1.000 

1.000 


1.000 


1.500 


1.500 


1.697 

1.697 

1.697 

1.697 

1.697 


1.697 


1.6 


1.6 


1.6 

1.6 

1.6 

1.6 

1.6 


1.6 


1.6 


1.6 


Voltage 

(V) 

17.5 

12.5 
12.5 


15.0 


12.5 


20.0 


Ibit 

(nN-s) 

107.00 

96.50 

79.00 


107.00 


96.00 


157.00 


106.0 

134.0 

96.0 

95.2 

73.0 


70.00 


66.20 


62.20 


61.1 

54.8 

52.3 

60.70 

68.20 


77.10 


67.20 


109.00 
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81. 

1.2 

12.5 

62.10 

43 

82. 

1.2 

15.0 

41.10 

43 

83. 

1.2 

17.5 

27.40 

43 

84. 

1.2 

20.0 

58.70 

43 

85. 

1.000 

1.2 

12.5 

58.00 

88 

86. 

1.000 

1.2 

15.0 

42.90 

88 

87. 

1.000 

1.2 

17.5 

32.00 

87 

88. 

1.2 

20.0 

57.40 

86 

89. 

1.500 

1.2 

12.5 

68.50 

130 

90. 

1.500 

1.2 

15.0 

38.70 

128 

91. 

1.500 

1.2 

17.5 

129 

92. 

1.500 

1.2 

20.0 

62.60 

131 

93. 

2.000 

1.2 

12.5 

70.20 

173 

94. 

2.000 

1.2 

15.0 

91.70 

175 

95. 

2.000 

1.2 

17.5 

47.00 

174 

96. 

2.000 

1.2 

20.0 

92.30 

175 
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Appendix  C  Thrust  vs.  Voltage  for  1.7  thru  1.2J 

This  first  graph  Figure  51  shows  how  the  thrust  varies  from  12.5  volts  to  20  volts,  the 
trend  is  that  there  is  a  slight  increase  in  the  thrust  levels. 


■  Thrust  at  different  voltages- 1. 47 J 


12.5  90.49762383 


0  50  100  150  200 

Thrust  (pN) 

Figure  5 1 .  Thrust  vs.  Voltage  for  1 .47J 


The  next  plot  is  the  1.697J  plot  and  it  has  different  behavior  and  not  the  smooth 
progression  that  is  was  in  Figure  51.  Here,  the  15  volt  setting  actually  has  the  highest  thrust 
number. 


i  Thrust  at  different  voltages- 1.697J 


> 


<D 

&Q 


12.5 

15 


£  17-5 


20 


70 


II  III  II  M  1  M 

72.47896707 

II  1  1  1 

72 


73.96445753 


73.61133687 

|  I  I  I  I 


76.93367449 


74 

Thrust  (pN) 


76 


78 


Figure  52.  Thrust  vs.  Voltage  for  1.697J 
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The  1 .6J  data  is  different  from  the  two  previously. 


i  Thrust  at  different  voltages- 1 ,6J 


> 


12.5 


15 


<D 

bD 

|  17.5 

> 


20 


30 


I  I  I  I  I 

138.14892806 


I  35.19081707 


32.91664538 


37.65180304 

I  I  I  I  '  I 


32 


38 


34  36 

Thrust  (pN) 

Figure  53.  Thrust  vs.  Voltage  for  1.6J 
i  Thrust  at  different  voltages- 1 .4J 


40 


65 


70 


85 


87.42890436 

I  I  I  I  | 

90 


75  80 

Thrust  (pN) 

Figure  54.  Thrust  vs.  Voltage  for  1.4J 

The  1.4J  data  appears  to  follow  the  same  path  as  the  1.47J  with  regard  to  increasing  thrust 
with  voltage. 
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The  thrusts  are  very  close  to  each  other  in  this  set,  probably  no  statistical  difference 
between  the  thrust  values. 


i  Thrust  at  different  voltages- 1. 3  J 


II  II  II 

oj.joJUjozy 

78.90184042 


20  40  60 

Thrust  (pN) 


80 


100 


i  Thrust  at  various  voltages- 1. 2  J 


I  I  I  I 

64.68945835 


20 


60 


67.74703737 

J-L-M 

80 


40 

Thrust  (pN) 

Figure  55.  Thrust  vs.  Voltage  for  1.3J  and  1.2 

The  thrust  appears  to  have  a  small  'bowl'  in  the  middle  of  the  voltage  ranges,  where  the 
voltages  are  less  than  the  lower  and  higher  end  values. 
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